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ABSTRACT 


Motivated by the study of genomes evolving by reversals, the primary topic of this 
thesis is “successful pressing sequences” in simple pseudo-graphs. Pressing sequences 
where first introduced by Hannenhali and Pevzner in 1999 where they showed that 
sorting signed permutation problem can be solved in polynomial time, therefore 
demonstrating that the length of a most parsimonious solution to the genome in- 
version only rearrangement problem can be determined efficiently. 

A signed permutation is an integer permutation where each entry is given a sign: 
plus or minus. A reversal in a signed permutation is the operation of reversing a 
subword and flipping the signs of the subword’s entries. The primary computational 
problem of sorting signed permutations by reversals is to find the minimum number 
of reversals needed to transform a signed permutation into the positive identity per- 
mutation. Hannenhalli and Pevzner showed that the signed sorting problem can be 
solved in polynomial-time in contrast to the problem of sorting unsigned permuta- 
tions, which is known to be NP-hard in general. At the core of the argument given 
by Hannenhali and Pevzner is the study of successful pressing sequences on vertex 
2-colored graphs. 

The connection between permutation sorting and phylogenetics dates back to at 
least the 1930’s, when two biologists, Dobzhansky and Sturtevant, wrote a series of 
papers in which they argued that the relationships between possible gene arrange- 
ments within a given chromosome encode critical information about the evolutionary 
history of species containing those genomes. In particular, they introduced the idea 


that the degree of disorder between the genes in two genomes is an indicator of the 


vi 


evolutionary distance between two organisms. This has inspired extensive work in 
the fields of computational biology, bio-informatics and phylogenetics. In particular, 
researchers have pursued the question of how a common ancestral genome may have 
been transformed by evolutionary events into distinct, yet homologous, genomes. In 
mathematics and computer science, we often represent genomes as signed permuta- 
tions (signed since DNA is oriented between two strands) and evolutionary events are 
encoded as operations on signed permutations. Among the most studied operations 
are block transpositions, prefix-reversals, and reversals, all of which correspond to 
common evolutionary mechanisms. 

In addition to the study of pressing sequences in simple pseudo-graphs, in this 
thesis we discuss related topics such as Cholesky factorizations of matrices over finite- 
fields, a sampling algorithm to generate simple pseudo-graphs uniformly at random, 
and the complexity of the “pressing space” of a simple pseudo-graph (the space of all 
successful pressing sequences of a simple pseudo-graph). This work includes collab- 
orative work with Dr. Joshua Cooper (Mathematics, University of South Carolina), 
M.S. graduate Erin Hanna (Mathematics, University of South Carolina), and M.S. 


candidate Peter Gartland (Mathematics, University of South Carolina). 
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CHAPTER 1 


ORIGIN OF PRESSING SEQUENCES IN GRAPHS 


1.1 INTRODUCTION 


In the late 1930’s, two biologists, Dobzhansky and Sturtevant, wrote a series of pa- 
pers in which they argued that the relationships between possible gene arrangements 
within a given chromosome encode critical information about the evolutionary history 
of species containing those genomes (see, e.g., [12, 27]). In particular, they introduced 
the idea that the degree of disorder between the genes in two genomes is an indica- 
tor of the evolutionary distance between two organisms. This has inspired extensive 
work in the fields of computational biology, bio-informatics and phylogenetics (see, 
e.g., [14]). In particular, researchers have pursued the question of how a common an- 
cestral genome may have been transformed by evolutionary events into distinct, yet 
homologous, genomes. Genomes are sets of chromosomes and chromosomes are built 
from deoxyribonucleic acid (DNA). DNA is a double-stranded molecule where each 
strand is oriented (corresponding to the direction that the ribosome reads the DNA) 
in opposite order. This means that in any of the DNA strand you have a footprint 
of all genes on the chromosomes: a gene is either in the order of the DNA strand 
(meaning the ribosome reads the gene from this strand) or the ribosome reads it from 
the other strand and what you see on the chromosome is the reverse-complement of 
the gene sequence. For this reason, in mathematics and computer science, we often 
represent genomes as signed permutations (the signs corresponding to the orienta- 


tions of the strands) and evolutionary events are encoded as operations on signed 


permutations. Among the most studied operations are block transpositions (e.g., [2, 
17, 30]), prefix-reversals (e.g., [7, 8, 15, 18]), and reversals (e.g., [1, 23, 16, 17, 20, 26, 


28]), all of which correspond to common evolutionary mechanisms. 


1.2. GENOME MUTATIONS 


The genome of a species can be thought of as a set of ordered sequences of genes. 
Each molecule is called a chromosome, and the set of all chromosomes is what we 
will call the genome. Chromosomes are made of deoxyribonucleic acid (DNA), a 
double-stranded molecule in which each strand is a long succession of nucleotides. 
DNA is made up of two complementary strands. During replication, these strands 
are separated and two copies of DNA are produced. Cellular proofreading and error- 
checking mechanisms ensure that the results are near-identical, but allow for some 
minor point mutations. There are three different kinds of point mutations: insertions, 
deletions, and substitutions. An insertion occurs when a nucleotide is inserted (added) 
to the sequence; a deletion occurs when a nucleotide is removed (subtracted) from 
the sequence; a substitution occurs when a nucleotide is removed from the sequence 
and replaced with another. 

A large-scale genetic mutation is referred to as a genetic rearrangement. While 
point mutations happen with some frequency during replication and typically have 
a minor effect on the outcome of the organism, genetic rearrangements happen in- 
frequently, have a large effect on the outcome of the reproduction, and therefore 
are a common metric used to estimate the evolutionary distance between organisms 
({12, 27]). The most commonly studied rearrangements include deletions, transposi- 
tions, inversions, duplications, reciprocal translocations, fusion, fission, and horizontal 


transfer. ([14]). 


e Deletions: A segment of the genome is lost. 


e Translocation: A segment of the genome moves to another location. 
e Inversion: A segment of the genome is reversed and the strands are exchanged. 
e Duplication: A segment of DNA is copied and inserted in the genome. 


e Reciprocal translocation: two broken off chromosomal pieces (each containing 


a telomere) are exchanged. 
e Fusion: Two chromosomes are joined into one. 
e Fission: One chromosome splits into two. 


e Horizontal transfer: A segment of the genome is copied from one genome to 


another. 


Originating from Greek words meaning same (homo) and proportion (logos), the 
term homology is use to describe the relationship between two objects that are similar 
in position or structure, but not necessarily in function. In mathematics, homology 
is a general way of associating a sequence of algebraic objects to other mathematical 
objects such as topological spaces. In genetics, two objects are said to be homol- 
ogous when their genetic sequences derive from a common origin. Observe that in 
the inversion model no genetic information is deleted or introduced, it is simply re- 
arranged. Thus, in this model, we consider two genomes to be homologous when 
they who present the same genetic information but in different orders. We attempt 
to understand how the two organisms may have evolved from a common ancestor by 
inversion rearrangements. 

In its simplest form, the genome rearrangement problem is formulated as follows: 
given a pair of genomes and a set of possible evolutionary events, find a most likely 
set of events that transforms one genome into the other. This thesis is motivated 
by studying the genome rearrangement problem restricted to inversions. Since rear- 


rangements are relatively rare events, scenarios minimizing their number are more 


likely to be close to reality ({14]). Thus, we may rephrase the genome rearrangement 
problem as such: given a pair of genomes and a set of possible evolutionary events, 
find a shortest set of events transforming one genome into the other; in particular, in 
the inversion model we are tasked with finding a shortest sequence of inversions that 
transform one genome into another. 

In the following sections we describe breakpoint graphs, and their corresponding 
overlap graphs. We will introduce successful pressing sequences on the bicolored 
graphs. In [20], the authors use successful pressing sequences on overlap graphs 
to demonstrate that the genome rearrangement problem restricted to inversions is 
solvable in polynomial-time. In particular, if we let N(A, B) be the minimum number 
of inversion operations need to transform genome A into homologous genome B, then 
for some integer k, N(A,B) can be computed in O(n") where n is the number of 
nucleotides in either genome. However, knowing N(A, B) does not necessarily tell us 
which sequence of N(A, B) edits can transform A into B; nor does it tell use how 
many such sequences exist. This question has been explored in [4] and in [9], and is 
the entry point for this thesis were we continue to explore this question as well as 


some related questions about sampling and complexity. 


1.3. BREAKPOINT GRAPHS 


In this section we introduce breakpoint graphs. This construction comes from [20]; a 
simplified version of the argument can be found in [3]. We begin by observing that, 
in the restricted view of evolution via reversals, a pair of genome can be considered 
homologous only if they contain the same set of nucleotides but appearing in different 
orders and in different orientations. Thus, given a pair of homologous genome with 
n nucleotides, we can label the nucleotides of one with the positively-signed identity 
permutation: 


id = (+1, 4+2,...,+n) 


and the other with a signed permutation 
Oo = (+01, +02, Sieg Eg) 


where o; is positive if and only if the nucleotide it represents is oriented in the same 


direction as the nucleotide that corresponds to it in the identity labeled genome. 


Example 1.1. Suppose we have a (very simplified) pair of homologous genomes, A 


and B, as illustrated in Figure 1.1. 


> 


Figure 1.1 Homologous Genomes 


We will label the 4 entries of genome B by (+1,+2,+3,+4). Genome A is then 
represented with the signed permutation (+4,—1,+2,+3). The genome rearrange- 
ment problem is to find a minimum number of evolutionary events (in this case 
reversals) to transform A into B, let that number be N(A, B). The following set of 


permutation reversals demonstrate that N(A,B) <5: 
(+4,-1,+2,+3) ~» (+4,-2,+1, +3) 
(+4,-2,+1,43) ~» (-1,4+2,—4,+3) 
(=1,+2,—-4,+3) ~» (+1,+2, —4, +3) 
(+1,+2,—-4,4+3) ~» (+1,+2, —4, -3) 
(elie 4, 8) (Gl 22 43,-4) 


The previous example illustrates that it is possible to upper-bound the number 


of inversions needed to transform one genome into another, however, these upper- 


bounds may be far from optimal. This lead to the creation of overlap graphs. We 
begin first by replacing each entry of o with a pair of entries according to the following 
rule: if 0; is positively signed then replace o; with 2-0; —1,2-0;, otherwise, when o; 


is negatively signed we replace o; with 2-0;,2-0;—1. 


Example 1.2. Suppose we have the same genomes as in the previous example. 
Genome A was represented with (+4, —1,+2,—3) is now represented with 

(7,8, 2,1,3,4,6,5). Genome B was represented with (+1,+2,+3,+4) is now repre- 
sented with (1,2,3,4,5,6,7,8). The reversal process in the previous example now 


becomes: 


(7,8, 2,1,3,4,5,6) ~» (7,8,4,3,1,2, 5,6) 
(7,8, 4,3,1,2,5,6) ~» (2,1,3,4,8,7,5,6) 
(2,1,3,4,8,7,5,6) ~» (1,2,3,4,8,7, 5,6) 
(1,2,3,4,8,7,5,6) ~» (1,2,3,4,8,7,6,5) 


(1,9; 3:48.76 ,5): ee (159) 3/45, 6,78) 


Observe that in every reversal the pairs of the form {2k — 1, 2k} are not separated at 
any point. The permitted reversals now are any reversals over a contiguous interval 


that has an even number of permutation entries to the left and to the right of it. 


We will now proceed to describe how to build the breakpoint graphs from a pair of 
permutations. Let o and 7 be unsigned permutations. By relabeling we may assume 


that 7 is the identity. 
me (Doo Oi): and.’ GOO 5, Cor) 


We will construct a plane graph (a graph embedded into the plane) with vertex set 
V = {0,1,2,...,2n,2n+1} drawn, in lex order, along a straight line (say the z-axis). 
Let A = {{2k,2k+1}|0<k <n}and B= {{oox, oon41} | O< k <n} where oo = 0 


and don41 = 2n+ 1. Let the edge set of the graph be EF = AU B, draw the edges 
of A as arcs above the z-axis with constant curvature, and draw the edges of B ina 


rectilinear fashion. 


Example 1.3. Consider the permutations o = (7,8, 2,1,3,4,5,6) and 7 = 

(1,2,3,4,5,6,7,8). The breakpoint graph of o and 7 (Figurel.2) has vertex set 
{0,1,...,9}, rectilinear edges {{0,7}, {8, 2}, {1,3}, {4,5}, {6,9}} and curved arcs 
{{0, 1}, {2,3}, {4, 5}, {6, 7}, {8,9}}. Observe that arcs aj = {2,2 +1} and ag = 
{y,y + 1} will cross in the breakpoint graph if exactly one of x and x + 1 appear 
between y and y+1 in o. Equivalently, if exactly one of y and y+1 appear between x 
and «+1 ino. Any proper drawing of the breakpoint graph will have 4 arc crossings 


(crossings of the interior portion of the arcs). 


Figure 1.2. The Breakpoint Graph of o and rT 


The problem of sorting signed permutations can be translated into performing a 
sequence of operations on the breakpoint graph which transform it into the identity 
breakpoint graph which has only trivial cycles. Before proceeding we introduce some 


terminology associated with breakpoint graphs. 


Definition 1.4. A breakpoint in a permutation occurs whenever two adjacent ele- 


ments are non-consecutive. 


Definition 1.5. The cycle number of a breakpoint graph, c(G), is the number of 
cycles in any proper drawing of the breakpoint graph. The breakpoint graph in 


Figure 1.2 has cycle number c = 2. 


Definition 1.6. The orientation of an arc is determined by the parity of the sum of 
their respective position. In particular, the arc a = {x,y} is positively oriented if the 
sum of the number of vertices drawn to the left of x and y is even; otherwise we say 


that the arc is negatively oriented. 


Observe that the goal of sorting signed permutations is to eliminate all breakpoints 
(implying that all the edges are positively oriented) by using the least number of 
reversals. In [2], Bafna and Pevzner showed that the reversal distance, d(G), in a 
breakpoint graph (the minimum number of reversals needed to transform the graph 
into the identity graph) satisfies the inequality d(G) > |V(G)|+1—c(G). Equivalently, 
d(G) > b(G) — ¢(G) where c'(G) is the number of cycles that contain four or more 


vertices. 


Definition 1.7. A cycle in a breakpoint graph is said to be an hurdle if all of its arcs 


are negatively oriented. 


In [20], Hannenhali and Pevzner showed that the reversal distance, d(G), in a 
breakpoint graph is 
d(G) = W(G) — G) + h(G) + f(G) 


where b(G) is the number of breakpoints in the non-identity permutation, c(G) is 
the cycle number of of the graph, h(G) is the number of hurdles in the graph, and 
f(G) € {0, 1} is an indicator of he number of fortresses in a graph (we will not define 
fortresses here except to say that they are a sort of hurdle that protects another 
hurdle). The argument that Hannenhali and Pevzner used relies on the construction 


of overlap graphs from breakpoint graphs. 


1.4 OVERLAP GRAPHS 


From a breakpoint graph we may construct an overlap graph in the following man- 
ner. Let B = (V(B), E(B)) be a breakpoint graph with V(B) = {0,1,...,2n+ 1}. 
Consider that B is already embedded into the plane. Let G = (V, EF) be a bicol- 
ored graph with V = {vo,v1,...,Un} and {u,;,u;} € FE if arcs a; = {2%,2¢ +1} and 


a; = {2j,2j7+1} cross in the drawing of B. Color v € V blue if degg(v) = 1 (mod 2). 


Example 1.8. Consider the permutations o = (7,8,2,1,3,4,5,6) and 
T = (1,2,3,4,5,6,7,8). Let B be the breakpoint graph of o and 7 (see Figure 1.2). 


The overlap graph of o and 7 is illustrated in Figure 1.3. 


M 


< 
Ss 


V, ON 


Figure 1.3. The Overlap Graph of 
o and Tt 


In the following chapter we will introduce simple pseudo-graphs which are the 
generalization of overlap graphs. The operation of pressing in a simple pseudo-graph 
corresponds to making a reversal in a breakpoint graph, and a successful pressing 
sequence corresponds to a most-parsimonious sequence of reversals that transform on 


permutation into the other. 


CHAPTER 2 


UNIQUELY PRESSABLE GRAPHS 


2.1 INTRODUCTION 


A signed permutation is an integer permutation where each entry is given a sign: 
plus or minus. A reversal in a signed permutation is the operation of reversing a 
subword and flipping the signs of the subword’s entries. The primary computational 
problem of sorting signed permutations by reversals is to find the minimum number 
of reversals needed to transform a signed permutation into the positive identity per- 
mutation. Hannenhalli and Pevzner showed that the signed sorting problem can be 
solved in polynomial time (see [3, 20]) in contrast to the problem of sorting unsigned 
permutations, which is known to be NP-hard in general (see, e.g., [6]). At the core 
of the analysis given in [20] is the study of “successful pressing sequences” on vertex 
2-colored graphs, which we refer to as “bicolored” graphs. In [9], the authors dis- 
cuss the existence of a number of nonisomorphic such graphs which have exactly one 


pressing sequence, the “uniquely pressables". In this paper we use a combination of 


graph theory and combinatorial matrix algebra over F2 to characterize and count the 


set of uniquely pressable bicolored graphs. 

This topic originated in computational phylogenetics, where Hannenhalli and 
Pevzner showed that certain bicolored graphs correspond to pairs of genomes and 
that the reversal distance between these genomes is the minimum length of a success- 
ful pressing sequence of said graph. In this interpretation, uniquely pressable graphs 


correspond to pairs of genomes linked by a unique minimum-distance putative evo- 


10 


lutionary history. The connection between permutation sorting and phylogenetics 
dates back to at least the 1930’s, when two biologists, Dobzhansky and Sturtevant, 
wrote a series of papers in which they argued that the relationships between possible 
gene arrangements within a given chromosome encode critical information about the 
evolutionary history of species containing those genomes (see, e.g., [12, 27]). In par- 
ticular, they introduced the idea that the degree of disorder between the genes in two 
genomes is an indicator of the evolutionary distance between two organisms. This 
has inspired extensive work in the fields of computational biology, bio-informatics and 
phylogenetics (see, e.g., [14]). In particular, researchers have pursued the question of 
how a common ancestral genome may have been transformed by evolutionary events 
into distinct, yet homologous, genomes. In mathematics and computer science, we 
often represent genomes as signed permutations (since DNA is double-stranded and 
hence oriented) and evolutionary events are encoded as operations on signed permu- 
tations. Among the most studied operations are block transpositions (e.g., [2, 17, 
30]), prefix-reversals (e.g., [7, 8, 15, 18]), and reversals (e.g., [1, 23, 16, 17, 20, 26, 
28]), all of which correspond to common evolutionary mechanisms. 

We continue the study of sorting by reversals by investigating its graph-theoretic 
generalization, “pressing sequences”, in graphs. Previous work in the area (see [4, 
9, 14, 20]) has employed the language of black-and-white vertex-colored graphs in 
discussing successful pressing sequences. For various reasons (such as simplifying 
definitions and notation), we find it more convenient to replace the black/white 
vertex-coloring with looped/loopless vertices. Thus, the object of study will be 
simple pseudo-graphs: graphs that admit loops but not multiple edges (sometimes 
known as “loopy graphs”). However, for the purposes of illustration we borrow 
the convention from [4, 9] and elsewhere that the loops of a simple pseudo-graph 
are drawn as black vertices. Given a simple pseudo-graph G, denote by V(G) 


the vertex set of G; E(G) C V(G) x V(G), symmetric as a relation, its edge set; 
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and G[S] = (S,(S x S)M E(G)) the induced subgraph of a set S C V(G). Let 
N(v) = Ne(v) = {w € V(G) : vw € E(G)} the neighborhood of v in V(G). Observe 


that v € N(v) iff v is a looped vertex. 


V; 


@-O O 


Figure 2.1. A simple pseudo-graph: 
G= ({v1, V2; us}, {v1v1, v1 03}) 


After this introduction, the discussion is arranged into four sections. In Section 


2, we develop some terminology and notation, and give a useful matrix factorization 


which we refer to as the “instructional Cholesky factorization” of a matrix (over F3), 
and discuss some of its properties. In Section 3, we present our main result, Theo- 
rem which characterizes the uniquely pressable graphs as those whose instructional 
Cholesky factorizations have a certain set of properties. In Section 4 we explore some 
consequences of the main theorem, such as the existence of a cubic-time algorithm for 
recognizing a uniquely pressable graph, a method for generating the uniquely press- 
able graphs by iteratively appending vertices to the beginning or end of a pressing 
sequence, and a counting argument which shows that there exist, up to isomorphism, 
exactly (3 — (—1)")/2-3!'"/2J-! uniquely pressable graphs on n non-isolated vertices. 
In the final section we discuss some open questions in this area. Before proceeding 
to Section 2 we list some basic notation for later use. Other terminology/notation 


employed below can be found in [5] or [11]. 


e We often write xy to represent the edge {x,y} for concision. In particular if 


x =y then zy is a loop. 
e [n] := {1,2,...,n} and [k,n] :={k,k+1,...,n} for all k,n EN. 
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e When S = V(G)\ {x} we write the induced subgraph of G on S, G[S], as G—2. 
In general, G — S' denotes G[V(G) \ SI. 


e For integers x and y, x = y (mod 2) is abbreviated as x = y. 


e For a square matrix M with rows and columns identically indexed by a set_X, 


for all x € X, Mz denotes the submatrix of M with row and column x removed. 


e When {a)}yca C Fo UZ, the notation >\¢, x, denotes addition over Z of Z), 


where Z), = x) if x, € Zand Z) is the least non-negative integer representation 


of x) in Z if x, € Fy. When referring to addition modulo 2 we use symbols © 


and @. For example, if 7] = 72 = 1 € Fy and x3 = 3 € Z then 


3 3 
Soa =1414+3=5 and Pej,=10101=1. 


i=1 i=1 


2.2. PRESSING AND CHOLESKY ROOTS 


Definition 2.1. Consider a simple pseudo-graph G with a looped vertex v € V(G). 
“Pressing vu” is the operation of transforming G into G’, a new simple pseudo-graph 


in which G|.N(v)] is complemented. That is, 
ViG)=V(G),  E(G!) = E(G@)A(N(v) x N(v)) 


We denote by Gi) the simple pseudo-graph resulting from pressing vertex v in V(G) 


and we abbreviate G(y,)(v9).--(o,) tO G(v,,09,...4,): For k > 1 we abbreviate (1,2,...,k) 


geeey 


as k so that when V(G) = [n] for some n > k then we may simplify Gi2,...4) to Gr. 
Go and Gy are interpreted to mean G. 
Given a simple pseudo-graph G, (v1, v2,...,v,;) is said to be a successful pressing 


sequence for G whenever the following conditions are met: 
e {v1,v2,--.,uR} CV(G), 


e v; is looped in Gy, v9... ) forall l <i<k, 


+Vi-1 
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. G(v1,v2,-.%e) = (V(G),9) 


In other words, looped vertices are pressed one at a time, with “success” meaning 
that the end result (when no looped vertices are left) is an empty graph. From 
the definition of “pressing” v we see that once a vertex is pressed it becomes iso- 
lated and cannot reappear in a valid pressing sequence. It was shown in [9] that if 
Gvr,02,.-.0n) = (V(G), 0) = Gor ut,...04,) then k = k’, ie., the length of all successful 
pressing sequences for G are the same. We refer to this length k as the pressing length 


of G. 


Definition 2.2. An ordered simple pseudo-graph, abbreviated OSP-graph, is a simple 
pseudo-graph with a total order on its vertices. In this paper, we will assume that 
the vertices of an OSP-graph are subsets of the positive integers under the usual 
ordering “<”. An OSP-graph G is said to be order-pressable if there exists some 
initial segment of V(G) that is a successful pressing sequence, that is, if it admits a 
successful pressing sequence (v1, V2,..., Ux) satisfying vu) < vg < +++ < vg and uz, < v' 
for all v' € V(G) \ {v1, v2,..., ve}. An OSP-graph G is said to be uniquely pressable 


if it is order-pressable and G has no other successful pressing sequence. 


Lemma 2.3. If G is a connected OSP-graph that is uniquely pressable then the press- 
ing length of G is |V(G)|. 


Proof. Suppose not. Then at some point in the pressing sequence two vertices become 
isolated and unlooped by one press. Call these vertices a and b, and suppose a is the 
vertex pressed. Then at this point in the pressing sequence a and b must have identical 


neighborhoods and hence 6 can replace a in the pressing sequence, contradicting that 


the pressing sequence is unique. 


We say a component of G is trivial if it is a loopless isolated vertex. 
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Proposition 2.4. /9/ A simple pseudo-graph G admits a successful pressing sequence 


if and only if every non-trivial component of G contains a looped vertex. 


Corollary 2.5. If G is a uniquely pressable OS'P-graph with at least one edge then G 


contains exactly one non-trivial component C and the pressing length of G is |V(C)|. 


Proof. Let G = ([n], £). Let C; and C, be (possibly distinct) non-trivial connected 
components of G. Let C3 = G—(V(C,) UV(C2)) so that G is the (possibly disjoint) 
union of C),C> and C3. As G is uniquely pressable it has unique pressing sequence 
og =n. Observe that pressing a vertex only makes changes to its closed neighborhood, 
a set which is contained within a single connected component. Let a; be the restriction 
of o to the vertices of Cj, 7 = 1,2,3. Then a; is a successful pressing sequence for 
G|[Ci]. If C: # C2 then G is the disjoint union of G[C)], G[C2] and G[C3], where the 
last one may be empty. Then pressing the vertices of G[C2] followed by pressing the 
vertices of G[C;] followed by pressing the vertices of G[C3] gives a successful pressing 


sequence for G’, contradicting the uniqueness of o. It follows that C; = C2 and o3 = 9, 


and therefore G contains exactly one non-trivial connected component. 


This shows that in order to understand uniquely pressable OSP-graphs, it suffices 


to understand connected, uniquely pressable OSP-graphs. 


Notation 2.6. CUP, is the set of connected, uniquely pressable ordered (<j) simple 


pseudo-graphs on n positive integer vertices. 


Definition 2.7. Given an OSP-graph G = ([n], £) define the adjacency matriz A = 


A(G) = (4,3) € F2*" by 


1 ifweebk, 
aij = 
0 otherwise. 


Note that A(G) is always symmetric. Previous work in the area refers to such ma- 


trices as augmented adjacency matrices as the diagonal entries are nonzero where 
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the vertices are colored black; since we have used looped vertices instead, the term 


“augmented” is not necessary. Define the instructional Cholesky root of G, denoted 


U= U(G) = (wes) — geen by 


1 ift<j andj € Na, _,(4), 
ij = ; 


0 otherwise. 

Observe that the j*” row of U is given by the j*” row of the adjacency matrix 
of Gj-1, and that u;, = 1 precisely when the act of pressing 7 during a successful 
pressing sequence of G flips the state of 7. Thus, U provides detailed “instructions” 
on how to carry out the actual pressing sequence. In Proposition 2.10 we justify use 


of the name Cholesky. 


Definition 2.8. For an order-pressable graph G = (|n], £) with instructional Cholesky 


root U = (u;,;) we define the dot product of two vertices i, 7 € V(G) as the dot product 


over F, of the i” and j*” columns of U: 


Ki 
(a; ja = B Ut iUt,j- 


t=1 


We define the (Hamming) weight of a verter j € [n] by 
wte(j) = 0 uy 
t=1 


and observe that wt¢(j) = (J.J) ¢- 


Figure 2.2. wt(1) = 1, wt(2) = 
wij 4, 12) 1, 3s) =), 
23) 1,24) =0,.3.4>=0 
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Definition 2.9. The weight of a column C in a matrix M, written wt);(C), is the 


sum of the entries in column C' (again, as elements of Z). 


Observe that if U is the instructional Cholesky root of G then the column weights 


of U correspond to the vertex weights of G. 


Proposition 2.10. Jf G = (|n], E) is an OSP-graph with successful pressing sequence 


1,2,...,k, adjacency matrix A, and instructional Cholesky root U then UTU = A. 


Proof. Let U'U = B = (b;,;) and A = (a;,;). Observe that b;; is the result (modulo 


2) of dotting the i” and j columns of U. Hence bj; = (i, j)a. 


For i,j € [n], let S;; = {t € [Kk]: 17 € E(Gr1z)AE(G;)} and T;,; = {t € [k]: t2 € 
E(Gt-1) and tj € E(G-1)}- 


Observe that S; ; lists the times during the pressing sequence that the pressed ver- 
tex results in the state of edge 77 being flipped. This occurs if and only if both 7 and 
j are in the neighborhood of the vertex being pressed. Hence S;; = T;,;. The state of 
the edge/non-edge ij in Gy is determined by its original state in G and by the number 
of times the state of the edge/non-edge 77 was flipped during the pressing sequence. 
However, G, = ([n],@) so ij ¢ E(G,) and therefore the number of times that the 
state of the edge/non-edge ij is flipped during the pressing sequence must agree in 
parity to with the original state of the edge/non-edge ij. It follows that |.5;;| = a;,;. 
On the other hand Tj, = {t € [n]: wi; = uz; = 1} list the common 1’s in columns 
i and j of the instructional Cholesky root. Hence |Jj,;| has the same parity as dot- 


ting the 7” and the j column of U. It follows that 6; = (7, 7)¢ = |Ti,j| = |Sij| = aij. 


Since the matrix entries are elements of F2, we have b;; = a;,; for 7,7 € [n] and 


therefore UTU = A. 


Le 


Observation 2.11. 
Given an OSP-graph G with adjacency matrix A and instructional Cholesky root U: 
ij € E(G) if and only if (i,7)¢ = 1, since a;,; is the result of dotting the i” and j™ 


columns of U. In particular 7 is looped in G if and only if wtg(7) = 1. 


In the theory of complex matrices, decompositions of the form A = U7U are 


known as “Cholesky” factorizations, so we repurpose this terminology here. While a 


symmetric full-rank matrix over Fy has a unique Cholesky decomposition (see [9}), 


a matrix M € F*" of less than full rank may have more than one Cholesky de- 


composition. On the other hand, the adjacency matrix A of an OSP-graph G with 
successful pressing sequence 1,2,..., has a unique instructional Cholesky root U as 
the first k rows are determined by the sequence of graphs G,G 1, Ge,...,G,_1 and 
the remaining rows (should they exist) are all zero. Throughout the paper we will 
take advantage of this by referring interchangeably to a pressable OSP-graph G, its 


(ordered) adjacency matrix A, and its instructional Cholesky root U. 


Example 2.12. Consider the F3*° matrix M and its (unique) instructional Cholesky 


root U: 
10 90620 4 10001 
010 10 01010 
M=j}0 00 0 0} and V=|0 000 0 
0. G0. 10 0000 0 
1000 1G 0000 0 


The following matrices also offer Cholesky factorizations for M: 


2D Oo Oo 
2D OF Oo 
So 6) e7 eo 
2D OF Oo 
FE Oo CO 
Sn 3S 3S 
2D OF Oo 
So 3S “Sr 7S 
2 OF Oo 
FE -F-& CO 
So +6 Oo eH 
oc OF oO 
So coco fo 
So OF oO 
oy He se: eS 
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Ke 
oS 
SS 
j=) 
me 
KS 
j=) 
j=) 
j=) 
mK 


Ot OL, 0 O20) 630 
000 1 0} and jo 0011 
000 1 0 00011 
000 0 0 0°00 070 


Definition 2.13. Consider a pressable OSP-graph G = (V, E) where V = {uj }ie{nj 
has the order implied by its indexing. For each j € [n] we say that v; has full weight 
in G provided wtg(v;) = j.. In particular if V(G) = [n] under the usual ordering then 
vertex j has full weight if and only if wtce(j) = 7 if and only if 7 € Ne,_, (2) for alli € 
[J]. 

The following notation will be used to simplify inductive arguments. 


Notation 2.14. For a given OSP-graph G = ((n], EF) with looped vertex j denote by 
GO = Gj) —j the result of pressing vertex j and then deleting it from the vertex set. 
Furthermore we let G? denote the result of pressing and deleting vertices 1,2,...,7 


in order from G. 


Lemma 2.15. Jf G € CUP,, has instructional Cholesky root U then Gt € CUP, _1 


and the instructional Cholesky root of G1 is U;. 


Proof. Let G = ({n],£) € CUP, with instructional Cholesky root U. The unique 


successful pressing sequence of G is n which is realized by 
G, Gi; Gay. rae (Gs 


and hence G; admits a successful pressing sequence: 2,3,...,n. Furthermore, if 
(v1, V2,---,Un—1) is a successful pressing sequence of Gy, then (1, v1, v2,...,Un—1) isa 
successful pressing sequence of G. By uniqueness it follows that v; = i+ 1 for each 
i €[n—1]. Then Gy admits exactly one successful pressing sequence 2,3,...,n, and 


therefore so does G!. It follows that G1 € CUP,,_;. Let V be the instructional 
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Cholesky root of Gt. The first row of G? is given by the neighborhood of 2 in Gy 


and in general the j*” row of V is given by N@1(23,...,3)(7 + 1). However 
Ner23,...3) (9 +1) = Na, (9 + 1) 


for each 7 € [n — 1]. Therefore the j row of V is the (j + 1)” row of U with the 


first entry deleted, since 1 is not a vertex in G!. V is the principal submatrix of U 


restricted to rows and columns 2,3,...,7. 


Proposition 2.16. /9/ An OSP-graph G = (|n],E) has pressing sequence n if and 


only if every leading principal minor of its adjacency matrix is nonzero. 


Lemma 2.17. Let G = (|n|, £) € CUP, with instructional Cholesky root U and let 
H = G-—n be the induced subgraph of G on |n— 1]. Then H € CUP,_; and the 


instructional Cholesky root of H is Up. 


Proof. If n = 1 then H = (0,0) which has only the empty sequence as a successful 
pressing sequence and its instructional Cholesky root is the empty matrix. Let n > 1. 
Observe that Ni(1) = Ne(1) — {n} so 1 is a looped vertex in H and therefore may 


be pressed to obtain Hy. For all 7 € [n — 1]: 


Na G)ANa(1), 19 € EA) 

Nu(J), 1j ¢ E(H) 

(NeG)ANe(1)) — tn}, ly € EU) @ 1j € E(G) 
Na(9) — (nh, 1j € EH) + 1j € E(G) 
Na,(j) —{n}, Ij € E(G) 


Na, (@j) — in}, 1 € EG) 
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Assume H; = G; —n for some 1 <i <n—1. Then 7+ 1 is looped in G;, implying 


that it is looped in Hj, so for all j € [n — 1]: 


Nu,G)ANa(i+1), {9,041} € EG) 
Naa ') 


(Ne (HANG, +1))— {nr}, {i,%+ 1} € E(Gi) 


Na,(9) — {nr}, {j,i +1} € E(Gi) 


Neu) = {n}, {3,4 a, 1} = E(G;) 


News) — {nr}, {t+ Ub € (Gs) 


By induction it follows that nm — 1 is a valid pressing sequence for H and H; = 


G; — {n} for all 7 € [n — 1]. We proceed to show that n — 1 is the only successful 
pressing sequence for H. Let A be the adjacency matrix of G (under the ordering 
n) and let U be its instructional Cholesky root. Let o = (v1, v2,...,Un—1) be a valid 
pressing sequence for H and let 7 = (v1, v2,.--,Un—1,n). Let P be the permutation 
matrix that encodes T. Then Aa is the adjacency matrix of H under the usual ordering 
< and P,A,P," is the adjacency matrix of H under the ordering given by a. Let V 
be the instructional Cholesky root of H under a. Observe that by Proposition 2.16, 
det(A) 4 0 and so det(PAP*) = det(P) det(A) det(P’) = det(A) 4 0. Furthermore, 


0 * 0 ‘ 
papt=| ms a 
0 * 0 
0 OF.) x x | 1 0 OF 4|,4 
* * 
P;,AaPa? VIV 
- * : * 
* * * * K * 
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Recall that H has n — 1 as a successful pressing sequence and so every successful 
pressing sequence must have length n — 1. It follows that all the diagonal entries of 
(upper/lower-triangular matrices) V and V* must be 1, implying that every leading 
principal minor of PAP® is non-zero. By Proposition 2.16, T is a successful pressing 


sequence for G. By uniqueness T = n and hence gd = n — 1. We may conclude that 


H € CUP,-1. 


Corollary 2.18. Let G € CUP, with instructional Cholesky root U. Then any prin- 
cipal submatrix of U on k consecutive rows and columns is the instructional Cholesky 


root of a CUP, graph. 


Proof. Follows by iteratively applying Lemmas 2.15 and 2.17. 


Corollary 2.19. If U is the instructional Cholesky root of GE CUP, then U must 


have all 1’s on the main diagonal and super-diagonal. 


Proof. Let H = ([2],£) € CUP,. Since it is connected, {1,2} € E; since it is 
order-pressable, 1 must be looped; and since it is uniquely pressable, Ny(1) #4 


Ny(2). Therefore CUP: = {([2], {{1, 1}, {1,2}})} which corresponds to instruc- 


1 1 
tional Cholesky root . The result holds by application of Corollary 2.18. 
0 1 


2.3 CHARACTERIZING UNIQUE PRESSABILITY 


Definition 2.20. For an upper-triangular matrix M € F3*” with columns C1, C2, 


..., C, with respective column weights w 1, w2,...,Wn, We say: 


cy = 1, We 


: Cj = (cay, C2jr+++ sy). has Property 1 if 
Cj =0, otherwise 


e M has Property 1 if each of its columns have Property 1. 


e M has Property 2 if 1 =w, < we <---< Wy 


Ze 


e M has Property 3 if w; > 2 implies wi42 > w;, for i € [n — 2]. 


e M has Property 4 if, whenever some non-initial column has odd weight, then it 


must have full weight and so must each column to its right. 


In other words, Property 1 is the condition that the nonzero entries in each column 
are consecutive and end at the diagonal; Property 2 is the condition that the weights 
of the columns are nondecreasing; and Property 3 is the condition that any column 
must have weight greater than that of the column two indices to its left if the latter 
has weight more than 2. Note also that Property 4 implies that any even-indexed 
column must have even weight; otherwise, it would have full weight, i.e., weight 
equal to the column index, which is even, a contradiction. Let M, = {M € F5*” | 


M satisfies Properties 1, 2, 3 and 4}. Observe that if Me M,, then Mp, € M,_1. 


Lemma 2.21. Letn > 1 and M € F5*” with columns and rows indexed by 1,2,...,n. 
If M € M,, then My € My-1. 


Proof. Let M = (¢i,3)ije{n] 80 that My = (ci; )o<ij<n- Let wi, We, ..,Wn be the column 


weights of 7. Let the columns and rows of Mj; be C2,...,Cn with weights w),...,w% 


n? 


respectively. Observe that w; =w,; — C1, for each 2 < 7 < n. It is immediate that 
M; inherits Property 1 from M. By Property 4 we know that the second column 
of M (as well as any even-indexed column of /) has even weight, it follows that 
W2 = 2 and so w) = 2 —c,2 = 1. Suppose towards a contradiction w; > wi,, for 


some 2<i<n-—1. Then 
wet] > wa +12 why tea = wa > wi > vj 
and so w; = w;. Then 
Wry < We = Wi S Wi = Wi +s S Wi +1 
and so wi41 = wi,, + 1. Hence we have 


Win = Wi, tl<wet+l=uw+1 > Witt S Uj. 
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It follows that w;,; = w; and therefore column 2 + 1 does not have full weight. By 
Property 1 of M we have c,;,; = 0 and therefore wi,, = wi4i, a contradiction. It 
follows that Mj has Property 2. 

Suppose now that 2 < wi for some 2 < 7 <n—2. Then w; > wi; > 280 Wj42 > W; 
by Property 3 of M. If w; = wj+2 — 1 then either column j or column j + 2 has odd 
weight which implies w;42 = j +2 by Property 4 of M. But then w; = 7 + 1 which 
is not possible. Therefore, 

WwW; Wael 


and 


which shows that M; has Property 3. To show Property 4 suppose w;, = 1 for some 


k > 2 (2 is the initial column of M;). Observe that 
w,t1> we > wy. 
If w, = w;, then w, = 1 and not full weight, contradicting Property 4. Hence 
Wy = Uw, +1 


and c;, = 1. It follows from Property 1 that w, = k > 2. Hence wi, > 1 and since 
w, = 1 then w), > 3 and wy, > 4. 

By applying Property 3 of M we get k+2 > were > we+1 =k +1. Since 
1=w, =k-—1 we have that wz. is the weight of an even-indexed column, and so 
Wee =O and wry =k 42. 


By arguing inductively, for all 7 € [|(n — k)/2]]: 


hence wy42; = k + 2). 


It follows that for all j € [|(n — k)/2]]: 
Wh49; = (k + 23) — C1 k427 = k + 23 — il. 
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Furthermore, for all 7 € [0,[(n—kK—1)/2]], Wriaj4 > Weta; = & + 27 and so 


Wr42j¢1 = k +27 + C1442541. Therefore, for all 7 € [0, [(n — k — 1)/2]]: 
Whey oj = Wepajt1 — Crktojp1 = k + 2). 


It follows that, in Mj, all the columns of index at least k have full weight and therefore 


M; has Property 4. 


Observe that the previous lemma can be extended to any matrix WM € M, by 
relabeling the rows and columns. We now proceed to our main theorem, which 
characterizes the set CUP,. This in turn provides a characterization of all the 
uniquely pressable simple pseudo-graphs (up to isomorphism), since the unique non- 
trivial, connected component of a simple pseudo-graph can always be relabeled to be 


a CUP graph. 
Notation 2.22. For an OSP-graph G let L(G) = {v | v € V(G) is a looped vertex}. 


Theorem 2.23. Let G = (|n|,£) with instructional Cholesky root U. Then G € 
CUP,, if and only ifU € My. 


Proof. For n = 1 the conditions of M, are only met by G = ({1], {(1,1)}) which in 
turn is the only full-length uniquely pressable OSP-graph on vertex set [1]. Let n > 1 
and assume towards an inductive argument that the statement holds for n — 1. We 
begin by showing sufficiency, that is if U € M, then G € CUP,. Choose and fix 
U = (uj) € M,. Let G = ([n], Z) be the OSP-graph with instructional Cholesky 
root U. By Properties 1 and 2, u;,; = 1 for each 7 € [n]. This implies that vertex 7 is 
looped in G;_1 for each 7 € [n]. It follows that n is a successful pressing sequence for 
G. We will show it is the only successful pressing sequence for G. Fix a successful 
pressing sequence o = 0j,...,0n for G. If o, = 1 then G(™ has adjacency matrix 


A(G!) = UTU; 
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By Lemma 2.21, U; € M,_; and therefore G() € CUP,,_; by the inductive 
hypothesis. Hence G‘“) has exactly one pressing sequence, and, since G(™) = G1, 
the sequence is (2,3,...,n). We may conclude that if 7, = 1 then o = n. Assume, by 
way of contradiction, that o, = t > 1. We will show that G contains a non-trivial 
loopless component and therefore is not pressable. Since t is a looped vertex it must 


have odd weight, and therefore full weight by Property 4. Let k = nin {i | wtg(i) = 
1} and let 


L=P,k= 1 RS lh) BSL). aid) Ra RUA 


By Property 4 of U, all the vertices in R have full weight. For all i € [n] and r € R: 


re = QP uniter = Bluar-1) 6 @B (us, - 0) = min {wtg(i), wtg(r)}. 


By Property 4, if wtg(?) = 1 then i € R. It follows that for all r € R, 


L(G) a[r—-1], ifr=o0 
L(G) U[r, nl, ifr=1 


Then 


However, (1,t) = 1 because t has full weight, so 1 € L(G) N Ne(t) and 
£( GM) CR: 


Similarly, for r € R 
New (r) = Ne(r)ANe(t) CR 


since 1 € Ne(r)N Ne(t). It follows that the induced subgraphs G'z) [| and Giz) {R| 
are contained not connected by a path in Gq. By applying Corollary 2.19, observe 
that 2 ¢ Ng(t) and 2 € Nc(1), so 


New (1) = Ne (IANc(t) > 2, 
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and therefore G) [z| contains an edge but no loops. It follows that GM does not 
admit a successful pressing sequence, a contraction. Therefore a, # 1. 

We now proceed to show necessity: if G € CUP,, then U € M,,. Let G € CUP, 
with instructional Cholesky root U = (u;,;). By Lemmas 2.15 and 2.17 we have that 
Gi1,G—n € CUP,,_; and therefore by the inductive hypothesis U;, U; € Mn_1. For 
simplicity we let H = G—n and w; = wtg(i) for i € V(G) throughout the rest of 


this proof. It suffices to show the following four conditions hold for U: 
(I) Property 1 holds for the n** column, 
(II). wy > Wat, 
(III) If w,z_2 > 2 then Wp > Wn_2, 
(IV) If w, #n then the first column of M is the only one with odd weight. 


We have four cases to consider. 


First Case: ui, = U2n = Ll. Since uz, = 1 then Property 4 of Uj gives us uj = 1 


for all 2 <i <n. It follows that w, =n and therefore (I), (II), (IID), and (IV) hold. 


Second Case: uj, = U2n = 0. (I) holds by Property 1 of U;. Recall that the 
diagonal and super-diagonal entries of U must be 1 by Corollary 2.19 so we need not 


consider the case where n < 3. For n = 4 we have two matrices to consider, 


1 10 0 1 1410 

OF 1. 20 Oe a tO) 
Y= and Va= . 

Oe OP at 001 1 

0001 0001 


V, satisfies (I) - (IV). Vo ¢ CUP, since (3,4,1,2) is also a successful pressing se- 


quence. Thus we may assume n > 5 to show (IJ), (III) and (IV) hold. 


yas 


Claim 2.24. u;,-1 = 0 


Proof of Claim 2.24. Assume towards a contradiction that ui,-1 = 1. Property 1 of 
U,, tells us that 


Wn-1 = wty, (n = 1) =n-1 


and so 


wty,(m — 1) = Wr. — Un = 1 — 2. 
By Property 4 of Uj, since uo, = 0, then 
wty,(n — 1) =0. 


It follows that 


Wn =n-1l=l. 


Recalling that ui, = U2n = 0, by Property 2 of U; we have 
n—2> wty,(n) > wty,(n—-1)=n—-2 


and so 


Let 


a min (2 | w; = 1}. 


Since G € CUP, and k 4 1 then G™) is not a pressable graph. We will use this to 


arrive at a contradiction. Since w, = 0, 


By the minimality of k, by Property 4 of Ua, and since n — 1 = 1: 


L(H) = {1} U {kk +2,...,n— 1}. 
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Observe that 


n 2 k n 
(kre = Butlin = @Buir- 0 CCAD) OS @ O- un) = k-2=1 
i=l i=1 i=3 i=k+1 


and for i € [n — 1] 
(k,t)g = (k, i)” = min(wty(k), wty(Z)) 
and so 
Na(k) = {1} U [k,n], 
It follows that 
L(G) = LIGJANg@(k) = {k 41,6 4+3,...,n—2,n}. 


Since N@(k)N|[k — 1] = {1}, the only potential edge in G [[k — 1]] affected by pressing 
k is the loop on 1. Furthermore G |[k — 1]] € CUP;_; by Corollary 2.18 so we may 


conclude that G |[k — 1]] is connected. If k A n—1 then (k+1,n)g =k—-1=0s0 
Ne(k +1) = {1,k} 
and so 
New (k +1) = Ne(k)ANe(k + 1) = [kK +1, n]. 
Then it follows that G™) [[k + 1,n]] is a connected graph with looped vertex (k +1). 
Observe that 
(l,n)a = 0, «l; ke = (k,n)a =1 

and so 


(1, N) Gtr) =] 


Therefore, if k 4 n — 1, G™ is a connected graph with at least one looped vertex, 
contradicting the fact that G € CUP,. We must conclude that k = n—1. Then 


Neo(k) = {1,2 —1,n} and so pressing k only affects four pairs of vertices: 
{(1,1), (1,7), (n —1,n —- 1), (n,n)}. 
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Once again, since G'[[k — 1]] is connected, G“) [[k — 1]] must be connected as well 


(although possibly without loops). However 
(n, n)e@ = 0, (1, ne = 0, and oe ka = (k,n)a =1 


sO 


CUnew = 1) and” “nyn)eey = 1: 


It follows that G is a connected graph with at least one looped vertex, namely n. 


This implies G is a pressable graph, contradicting that G € CUP,. Claim 2.24 is 


established. 


Claim 2.25. U1n-2 = 0 


Proof of Claim 2.25. Assume towards a contradiction that uj,-2 = 1. By Property 
1 of Us, 


Wn_2 = wty,(n — 2) =n—-2. 


By Claim 1 


wty,(n —1)<n-1 


and so by Property 4 


wty,(n — 2) =0 


and therefore n = 0. We then have 


wty,(n — 2) = Wn-2 — Uin-2 =N-3=1 


which contradicts Property 4 of U; since u2,, = 0. This establishes Claim 2.25. 


We may now assume Uj,-2 = Uin—1 = 0 and proceed to show (II)-(IV). (I) 


follows from Property 2 of Uj; since 


Wn = Wty, (n) = wty,(n — 1) = Wy-1. 
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To verify (III), observe that if w,_2 > 2 then 
wty,(n-2)>2 > wr=wty,(n) > wty,(n — 2) = wre 
by Property 3 on U;. (IV) is established by observing that 
Un1=0 => wty,(n-1<n-1 = wty,(i) =O forall2<i<n 
by Property 4 of U,; and so 
w; = wty, (i) =0 for all2 <i<n. 


Third Case: u;,, = 0, v2, = 1. 

By Property 1 of U; we have that u;,, = 1 for all2 <i <n. It follows that w, = n-1 
and so (I) and (II) hold. Furthermore, since w,_2 < n — 2 then (III) holds. To verify 
(IV), we need to show that L(G) = {1}. 


Claim 2.26. L(G) ¥ {1, n}. 


Proof of Claim 2.26. Assume, by way of contradiction, that £(G) = {1,n}. Since 
G € CUP, and n 4 1 then G™ must contain a non-trivial loopless component C. 
Choose and fix p € V(C). Since C is non-trivial we may assume p 4 1. Since H € 
CUP,,_: there must exist a path from p to a looped vertex in H. Since £(H) = {1} 


this looped vertex must be 1. Choose such a path P, 
P=v9...V¢ 
where p = Up and vg = 1. Observe that 
Clea a0 


so 1 is a looped vertex in G( as well. Then some interior edge of P must be removed 


upon pressing n as otherwise p would have a path to a looped vertex in G. Let 


j= mint? | (Me = 1p 


dl 


then 


P’ = uP; 


is a path from p to a looped vertex in G“?, a contradiction. This establishes Claim 


2.26. 


Claim 2.27. Ifn =0 then L(G) = {1}. 


Proof of Claim 2.27. Let n = 0. Assume, by way of contradiction, that L(G) 4 {1} 
and let 


k= juin | w; = 1}. 


Choose a non-trivial loopless component C of G*) and a vertex p € V(C) \ {1}. 


Recall that ui, = 0 and u2, = 1. Then 


(n,ne= Pun =n-1= 1 and (k,n)g = Bui, =k-1 = 0 


i=l i=l 
implies 
(n, N) otk) ale 


son #p. For any j € [k+1,n—1]\L(G), since H € CUP, _1 and w; = 1, Property 


3 implies that (k, 7) = 1, whence 
GIG = 0 and (k, J)e@ =1 implies (9,9) GH =1 


so j is looped in G™) and therefore p ¢ [k+1,n—1]\L(G). Ifj € [k+1,n—1]NL(G) 


then j € [k+2,n —1] and 
J = leg = 0 and (ky J = Ye = (kj )a@ =1 implies See) = Lew =1 
so p ¢ [k+1,n—1]NL(G). Therefore 


p€ [2,k-1). 
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Since G [[k — 1]] € CUP,_,; then there exists a path P in G'[[k — 1]] connecting p to 
a looped vertex; since £(G|[kA — 1]]) = £(G) n [k — 1] = {1} then the looped vertex 
must 1. 


PS 052202 


where up = p and ve = 1. Note that 
Ne(k) OV(P) = {1} 


because, if i = min{j : v; € Ng(k)} is not 1, then v; is looped in G™ and in the same 
component (namely, C’) with p, a contradiction. Then all the interior edges of P are 


unaffected by pressing k (although the loop on 1 is removed). If 2 € V(P) then 
P! = upP2n 

is a path from p to n in G“). If 2 ¢ V(P) then 
P’ = vpPy2n 


is a path from p to n in G). Since n is looped in G this contradicts that C is a 


non-trivial loopless component. Claim 2.27 is established. 


By Claims 2.26 and 2.27 we have only one case left to consider. Let n = 1 and 


assume, by way of contradiction, that L(G) 4 {1}. Let 
k= min AU | w; = 1}. 


31S n— 


Choose a non-trivial component C’ of G“) and a vertex p € V(C) \ {1}. Observe that 


n—1=0s0 wtg(n — 1) = 0 by Property 4 of H. Then, by Property 4 of H 
(n-1n-l)g=n-1=0and (kn-lIbg=k=1 


which implies 


(n-—1,n-l)gem =1 
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son—1+#p. Similarly, 


n-1 k 
(n—1,n)g=1-08@1-1l=n-2=1and (k,n)c =1-06@P1-1=k-1=0 
i=2 1=2 


implies 
(n-—1,n)em =1 
so n # p as it is connected to a looped vertex. If 7 € [k + 1,n — 2] \ L(G) then 
9)G = 0 and (k, ja =1 implies J) G@) =1 

sop¢ |kK+1,n—2]\ L(G). If7 € [k+1,n—2]N L(G) then j € [A +2,n — 2] and 

Ohi _ le = 0 and (ky j aa lhe = (k, j)e@ = 1 implies Ore) a lew =1 
so p ¢ [k+1,n—2]N L(G). It follows that 

p€ [2,k-1). 


Since G |[k — 1]] € CUP,_1 has only 1 as a looped vertex then G [[k — 1]] contains a 
path 


PH Vged Op 


where vp = p and vy = 1. Since Ne(k) N [k — 1] = {1} none of the interior edges are 


removed upon pressing k. If 2 € V(P) we have 
P’ = v9 P2n(n — 1) 
is a path to a looped vertex in G“) and if 2 ¢ V(P) then 
P’ = vpPwy%2n(n — 1) 


is a path to a looped vertex in G“). This contradicts that C is a loopless component 


in G“), We conclude that £(G) = {1} and therefore establish (IV). 
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Fourth Case: u;,,, = 1, U2, = 0. 
We show that this case cannot occur. Assume it does to reach a contradiction. Since 
Un = 0 then by Property 4 only the initial column of Uj has odd weight. It follows 
that 

Wn = Wty,(n) +121. 
Let k = gain (3 | w; = 1}. 
Claim 2.28. k#n 
Proof of Claim 2.28. Assume, by way of contradiction, that k = n. Then L(G) = 
{1,n}. Since G € CUP,, we may conclude that G contains a non-trivial loopless 


component C. Choose and fix p € V(C) \ {1}. Since H = G[[n — 1]] € CUP, -1 


there exists a path P in H that connects p to a looped vertex, namely 1. Let 
P=v9...V¢ 
where vp = p and vg = 1. Observe that 
onne=(Lne=1 


and let 


i= gaint: | vu; € Ne(n)}. 


If m < @ then 


P! = Pum 


is a path to a looped vertex in G), contrary to assumption. Assume m = @ and let 
q = Ue-1- 


tS (ve, Ve-1)G = (lade 


SO Uiq = 1. By Property 1 of Uj, wig = 1 for all 7 € [gq] and so wy = qg. Let 


c= gain (3 li ageacdips 
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Since (q,n)g = (vp_1,)g = 0 then 


0= Dlg) = (le I) Bis 0) @ ou -1)@ @ (0-1) 


and therefore 


q 
(n)a=0=10Q@1 


i 


which implies r < q. Furthermore g 4 n — 1 since otherwise we would have 
Wn-1 = Wg =GQ=an-—-1 
and 
n—-1> Un = Un + Wty, (n) = 1+ wty,(n) = 1+wty,(n—1) =1+(n—-2) 


which would imply 

Un =n-l=w,=0 
contradicting that n is a looped vertex. Thus q 4 n—1 and so g+1 <n which gives 
us that wg41 = 0 by Property 4 of U; and the fact that uz, = 0. However, 


g=w, = wty,(q) < wty,(¢q+)D) =wai <¢qtl 


implies wgi1 = gq since g+ 1= 1 and therefore u1,,4; = 0 by Property 1 of U;. Then 


n qd n 


(94+ No =Q@ig: vig) = (1-0 6G1-1)80-1)e @ (-0)=1 
i=l i=2 i=qt+2 

and 

r-1 gt+l n 
(q ne 1,n)e = (0 : es D Dui gt . 0) D Bi . 1) @ DB (0 tea) 

i=2 i=r i=qt2 

and so 
qtl q 


Observe that 


2 
(2, G+ l)c = BD wigs =1 


i=1 
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sO 
P'=q(q4+1)2 


is a path in G whose q(q + 1) and (q + 1)2 edges are unaffected by pressing n (since 


(¢+1) € Ne(n)) and since 
2 2 
(2,2)¢=@QMu2=0 and (2,n)\¢=Qun=1 


i=1 i=1 


P* = pPqP'2 
is a path from p to a looped vertex in G™, contrary to assumption that p is contained 


in a loopless component of G(™. This contradicts that G € CUP, and establishes 


Claim 2.28. 


We proceed under the assumption that 
k= gain (3 |w=lb<n, 


once again in search of a contradiction. Observe that we need not consider the case 
where n < 3 since the super-diagonal entries must be all 1. Furthermore, if n = 4 
we have only two matrices to consider, both of which have an additional successful 


pressing sequence given by (4, 3, 2, 1): 


1 101 1141éi1 
01410 01410 
and 
0011 0011 
000 1 000 1 


Assume n > 5. By Property 4 of Uj; we have u1,,-1 = 1 and so by Property 1 and 2 
of Us 


n—-1> Wp, = wty,(n) +1 > wty,(n-1) +1 =(n-2)+1=n-1 


which gives us w, =n —1. Observing that wty,(n) =n —2 and uz, = 0 we conclude 


from Property 4 of U; that n — 2 = 0, so n is looped in G. Since G € CUP, then 
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Gm) must contain a non-trivial loopless component, say C’. Choose and fix a vertex 
p € V(C) \ {1}. Observe that G [[p]] € CUP, by Corollary 2.18. If p € L(G) then 


pressing n must remove its loop so 
(2,p)a =0 and (2,n)¢ = (p, na =1 


which implies that (2,p)¢,,,) = land 2 € £(G(n)), contradicting that C is loopless. It 
follows that p ¢ £(G) and therefore, in G [|p|], we have a path P from p to a looped 
vertex b: 


P= Vpn icVp 


where v9 = p, vg = 6, and v; is loopless in G[[p]] for 0 < i < @ Observe that 
Ne(n) 9 {v1, v2,..., ve} = @ since otherwise we would have a looped vertex in C. It 
follows that the interior edges of P are unaffected by pressing n in G and therefore 
it must be the case that pressing n in G removes the loop from b = ve. By Property 


4 on Us looped vertices in H have full weight. Observe that 


Dui ptin) = (1-1) @ (24-0) € Buin: 1) = 
= = 1+(b-2), ifbf1 


Since 1 = (b,n)g = BD (wi pttin) this implies that b = 1. (Otherwise, b is even but 
i=l 


looped in G, contradicting Property 4.) Observe that wz; = 3 would imply w4 = 
wty, (4) = 4 by Property 4 of U,, and then wty,(4) = 1 but uz, = 0, contradicting 
Property 4 of U;. Then 


ts = wip, (3) = 2 
by Property 2 of U,; and so 
@B,ae¢=0 nd “Ges anplies:- <3,3) 6621 
and so 3 ¢ V(P). Furthermore 
1,3)¢=0 and (1,n)g=(8,n)g=1 implies (1,3)e@m =1. 
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Therefore 


is a path to a looped vertex in G™. This implies G ¢ CUP, contrary to assumption. 


Therefore Case 4 cannot occur. 


P’ = vo Pv,3 


2.4 RECOGNITION AND ENUMERATION 


A straightforward and very slow way to check if a simple pseudo-graph on n vertices 
is uniquely pressable is to check the pressability of each one of its n! orderings. Here 
we offer a substantially faster algorithm. 
Algorithm 1: Find a Pressing Order 

1: input: Adjacency matrix A with entries a;,; for i,j € [n] 


2: output: Re-indexed matrix P? AP 


3 McA 

AE t= Oasen 

ae Ae el 

6: while t < n do 

7: maxDegree + 0 

8: indexMax Degree + 0 
G. -4 eed 


10:  whilei<ndo 


11: deg; < 0 

12: if m;,;, = 1 then 

13: deg: <— Vij Mi,j 

14: if deg; > max Degree then 
1b: max Degree < d; 

16: index Max Degree < i 
17: Ue cee a 
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18: if indexMaxDegree = 0 then 


19: if Ve Vj-1 mej > 0 then 

20: return False {Not a Pressable Graph} 
21: else 

22 k + index Max Degree 

23: (Oe ae 

24: hah ae 

25: for ¢ € [n] do 

26: for 7 € [n] \ {k} do 

Pa Me; — Me; B (Mn; > Mex) 
28: for j € |[n] do 

29: ting =H 0 


30: return P! AP 


Algorithm 2: Construct instructional Cholesky root 
1: input: Adjacency matrix A with entries a;,; for i,j € [n] 
2: output: Instructional Cholesky matrix U 
go Sy 
Ake 1 
5: while k <n do 


6: if Ak,k = 1 then 


7: for j € |[n| do 

8: Uk,j *— Gk 

9: for i € [k+1,n] do 

10: for j € {[n| do 

11: Ai — Giz B (Akz - Gin) 
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12: ke k+l 
13: else 
14: kene+l1 


15: return U 


Algorithm 3: Does this instructional Cholesky correspond to a uniquely pressable 
OSP? 

1: input: Instructional Cholesky matrix U with entries u;; for 7,7 € [n] 

2: output: True or False 

3.9 1 

4: while j <n do 

Bo - 2421 


6: while: <7 do 


ce if uj —0 then 

8: (4-4-1 

9: else if J_, Uej <j —t+1 then 
10: return False 


ll: gegtl 

IR i ete 

13: while 7 < n do 

14: if 3. ues > WF) wey then 

ies return False 

16: else if j <n—2 and Y3_, ue; > 2 and D344 ue j4e = Di, wey then 
ine return False 

is: else if Jt} ue jar = 1 and D934) wey AJ +1 then 

19: return False 


20: else 


Al 


21: fd 


22: return ‘True 
Corollary 2.29. The unique pressability of G can be decided in time O(n). 


Proof. Let G = (V,E) be a simple pseudo-graph on n vertices. Let G’ = ((n], E’) 
be the result of relabeling of V(G) so that G’ is order-pressable. If G is uniquely 
pressable then the instructional Cholesky root of G’ is in M,,. Observe that for each 
ke L(G) \ {1}, 

Ner(k) = {1} U [k,n] 
since by Property 4 each @ € [k,n] has full weight and therefore (k,0)q@ = k = 1. 
However, for each k € L(G") \ {1} and @ € [k,n], (1, 2)q@ = 1 by Property 4, and 


ae 2)Gr =] and en k) ou =2=0. 


Thus, Nea (1) ») {2} U Na (k) 2 No (k). 
Therefore 1 is the unique looped vertex of largest degree. It follows that to 
find a (potentially unique) pressing order it suffices to iterate the process of finding 


the looped vertex of largest degree and pressing it. Index the vertices of graph G 


arbitrarily and define A = (a;,;) € F5*", the adjacency matrix of the graph. We 


proceed to describe three algorithms which have been included in the appendix. 
Algorithm 1 finds a successful pressing sequence for G (given that one exists) 
by finding the looped vertex of largest degree and pressing it; this has running time 
O(n°), as it amounts to performing in-place Gaussian elimination on an n Xn matrix. 
Algorithm 2 computes the instructional Cholesky root of an ordered adjacency ma- 
trix and once again is done by performing Gaussian elimination. Finally, Algorithm 
3 checks if an upper-triangular matrix has the properties of M,, which is done by 
computing no more than n partial sums for each of n columns and comparing them 


sequentially. Algorithm 3 also has running time O(n?). 
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Corollary 2.30. Let n > 0. Let G = ((n|, E) and let H = ([n + 1], E(A)) be the 
result of adding a vertex “n +1” adjacent to each looped vertex in G, with a loop at 


n+1 if and only ifn ts even. IfG € CUP,, then H € CUP, 41. 


Proof. Suppose G € CUP, with adjacency matrix A and instructional Cholesky root 
U = (uj). Observe that H = ([n + 1], E(H)) where 


E(G)U{(i,n +1) lie L(@}U{(nt1n+D}, ifn=o 


E(H) = 
E(G)U {(i,n +1) |i e L(G)}, ifn=l 
1 
U ‘ 
Let V = and observe that 

1 
0 > Q | 1 

Di n+1 

A 
YVvS 
Uae 
Daca eee | peas On+1n+1 


where bp41; = 1 if and only if wty(¢) = 1 if and only if7 € L(G) ort =n4+1= 1. 
It follows that V7V is a Cholesky factorization for the adjacency matrix of H. Since 
G € CUP, then U has 1’s along the diagonal and so U and V are full rank matrices. 
Since Cholesky factorizations are unique for full-rank matrices [9] then V must be 
the instructional Cholesky root of H. V inherits the properties of M,, in its first 


n columns from U. Furthermore the last column of V has full weight n + 1, so 


V € Mnii, and therefore H € CUP,41. 
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Corollary 2.31. Letn > 0. Let G = ([2,n+1],E£) and H = ([n+ 1], E(A)) be 
the result of removing all the edges, including loops, from G[L(G)|, adding a looped 
verter “1” adjacent to all of L(G). If Ge CUP, then H € CUP) 41. 


Proof. Let G € CUP, recall that by Property 4 (of its instructional Cholesky root) 
the looped vertices in G form a clique. The only looped vertex in H is 1 so if H 
admits a successful pressing sequence it must begin with 1. However H‘) = G since 


pressing and deleting 1 creates an edge between any two vertices in L(G). It follows 


that H has exactly one successful pressing sequence: n+ 1. H € CUP,41. 


Notation 2.32. CUP), is the set of connected, uniquely pressable ordered <y simple 


pseudo-graphs on vertex set [7]. 


Corollary 2.33. Up to isomorphism, the number of connected, uniquely pressable 


simple pseudo-graphs on n > 1 vertices is 


cup | Bh). n is even 
: 2. 3(-8)/2 nn Gs odd. 


Proof. For n = 2, the result holds since 


We proceed by induction. Let n > 2 be even and assume |CUP in = 3("-2)/2. Choose 
and fix G = ({n], £) € CUP), with adjacency matrix A and instructional Cholesky 
root U = (u;,;). Let G’ = ([2,n+ 1], E’) be a re-indexing of G given by it)» i+ 1 for 
all i € [n]. Let Hy = ([n + 1], £), Ho = ([n + 1], Ee) where 


BE, = EVU{(i,n+1) [ie L}u{(n¢1,n+)} 


and 


Ey = B'A ((£(G') U {1}) x (L(G) U {1})). 
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By Corollaries 2.30 and 2.31, Hy, Hz € CUP {n41);. Hy has at least two looped vertices 
(1 and n+ 1) and Hy has only one looped vertex (1), so H,; #4 Hy. Furthermore, the 
instructional Cholesky roots of H, and H» include U as a principal submatrix on 
consecutive rows and columns so it is not possible that we would have gotten H, or 


H by applying Corollaries 2.30 or 2.31 to other graphs in CUP,,). It follows that 
|CUP in1)| >2- |CUP in] = 9 BAP. 


Consider now any H € CUP ),+1; with instructional Cholesky root V = (u,,;). If 
Vint = 0 then let G’ = H™ and let G be the re-indexing of V(G) given by its i-1 
for alli € [2,n+1]. By Lemma 2.15, G € CUP),) and therefore H can be constructed 
from G using Corollary 2.31. If vin41 = 1 then viny1 = 1 for all 7 € [n + 1] because 
of Property 4 by appeal to Theorem 2.23. Let G = H — {n+ 1}. By Lemma 2.17, 
G € CUP, and hence H can be obtained by applying Corollary 2.30 to G. It follows 
that 
|CUP jn41]| < 2-|CUP pn] = 2-3-9”. 


We count |CUPin+2]| in a similar, though slightly more complicated, manner. Given 
G € CUP,,), let H, be result of two successive applications of Corollary 2.30. Let 
Hy be the result of applying Corollary 2.30 followed by Corollary 2.31, let H3 be 
the result of applying Corollary 2.31 followed by Corollary 2.30, and let H4 be the 
result of applying Corollary 2.31 twice successively. We first show that Hz = H3 and 
then argue that H,, Hy, H,4 are pairwise distinct. Let V2 and V3 be the instructional 


Cholesky roots of Hz and H3, respectively. Then 


1 bio ieee bin 1 
0 1 
Y= U = V3 
0 1 
0 QO. 0 1 
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where b; ; = 1 whenever it is positioned above a column of odd weight in U. Therefore 
Hy = Hz. Observe that H, has at least two looped vertices, 1 and n+ 1, whereas Hy 
and Hy, have only one looped vertex “1”. Since n is even, vertex n in G is loopless by 


Property 4. Then the instructional Cholesky root of H4 has the form 


1] 1 big -+* Diss 
0 if bo3 Aner bo n+2 
Ya= 0 |o 
U 
0 Jo 


where 01 n42 = ben42 = 0 and so Hy ¥ Ap. It follows that 


|CUP in| >3- |CUP i, — atar2)-3)/2) 


Choose and fix H € CUP),+2) with instructional Cholesky root V = (u,;). Let 
G, = (H —{n+2}) —{n+1}, G, = HY — {n +2}, and G5 = H%?). Let Gp and 
Gs be (order-preserving) re-indexings of G, and G5 so that V(G2) = V(Gs) = [nl]. 
By (repeated) applications of Lemmas 2.15 and 2.17; G,G2,G3 € CUP). Let 
=U ntit Vint Ifa = 2 then v7.41 = Viny2 = 1 and by Property 1, vin41 = 1 
for all i € [2 + 1] and U;n42 = 1 for all i € [n + 2]. Then A can be constructed from 
two applications of Corollary 2.30 to Gy. Ifa = 1 then v,4; = 0 and v4,42 = 1 by 
Property 4. Furthermore, Property 4 implies that 1 is the only looped vertex, since 
n+ 2 is even and v,,41 = 0. Observing that n+ 2 must be a full weight vertex, we 
conclude that H can be constructed from G2 by application of Corollaries 2.30 and 
2.31 (in either order). Finally if a = 0 then v1,,42 = 0 and by Property 4, and since n 
is even, V2,n42 = 0. Furthermore by Property 4 it follows that H and H™) have each 
only one looped vertex. Then H can be constructed from G3 by two applications of 


Corollary 2.31. It follows that 
ICUP ja42}| <3 -|CUP iy] = 34-9? 
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Therefore 


Bee Ae, if n is even 


|CUP)n| = . 
2-3-3)/2 if n is odd 


Corollary 2.34. The number of uniquely pressable simple pseudo-graphs onn > 1 
vertices up to isomorphism is 


(5 - gln—2)/2 4 1) /2, ifn is even 
7 = 


re? ~ 1) 12s if n is odd 

Proof. There are three non-isomorphic uniquely pressable simple pseudo-graphs on 2 
vertices: the edgeless (loopless) graph, the disconnected graph containing one looped 
vertex and one unlooped vertex, and the connected graph containing one looped 
vertex and one unlooped vertex. We proceed by induction on n. Observe that for 
every k <n and for every H € CUP,, we can create a (distinct) uniquely pressable 
graph G on n vertices by adding n — k isolated vertices to H. Similarly, if G is a 
uniquely pressable graph, then it is either the edgeless (loopless) graph or it contains 
exactly one non-trivial component which must be a CUP, graph for some k < n. 
Hence 


T, = ))|CUP,y| = Tr-1 + [CUP]. 
k=0 
The result follows by observing that 


n—1)—2)/2 n—3)/2 n—3)/2 n+1)/2 
5 . 3((r-1) P+l  gin-ay2_ 5:3! V2 44 44. 3(n-3)/ 3 W244 
2 2 2 


and 


3((n—D+1)/2 4 4 3. 3(r-2)/2 4 1 4.9. 3(n-2)/2 5s 3(n—2)/2 4 1 


1 3(0-2)/2 — ; ; 


AT 


CHAPTER 3 
AUTONOMOUS POSETS 


3.1 INTRODUCTION 


A simple pseudo-graph is a graph that admits loops but not multiple edges (sometimes 
known as a “loopy graph”). Given a simple pseudo-graph G, denote by V(G) the 
vertex set of G; E(G) C V(G) x V(G), symmetric as a relation, its edge set. Let 
N(v) = Ne(v) = {w € V(G) : vw € E(G)} the neighborhood of v in V(G). Observe 
that v € N(v) iff v is a looped vertex. For S C V, we denote by G[S] the vertex- 


induced subgraph on S. 


Definition 3.1. Consider a simple pseudo-graph G with a looped vertex v € V(G). 
“Pressing uv” is the operation of transforming G into G’, a new simple pseudo-graph 


in which G|N(v)] is complemented. That is, 
ViG)=V(G), E(G'!) = E(G)A(N(v) x N(v)) 


We denote by G'(,) the simple pseudo-graph resulting from pressing vertex v in V(G) 
and we abbreviate G(y,)(v2).-(vg) tO G(w,,v9,...,ve): For k > 1 we abbreviate (1,2,...,k) 
as k so that when V(G) = [n] for some n > k then we may simplify Gi2,..4) to 
Gy. Go and Gi are interpreted to mean G. To aid with inductive arguments, we let 
GG) = Gv) — uv: the result of pressing v in G (which leaves it isolated, loopless, and 


thenceforth unpressable) and then removing the pressed vertex. 


Given a simple pseudo-graph G,, (vj, v2,...,v;) is said to be a successful pressing 


sequence for G whenever the following conditions are met: 
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° {U1 W354 54 VEY c V(G), 
e v; is looped in Gm, 29,...,v,-1) for all 1 <2 < k, 
e Gai cas = (V(G),0) 


In other words, looped vertices are pressed one at a time, with “success” meaning that 
the end result (when no looped vertices are left) is an empty graph. This topic orig- 
inated in computational phylogenetics, where Hannenhalli and Pevzner showed that 
certain simple pseudo-graphs correspond to pairs of genomes and that the reversal 
edit distance between these genomes is the minimum length of a successful pressing 
sequence of said graph [20]. In phylogenetics, the simple pseudo-graph corresponds 
to a pair of homologous genomes and its successful pressing sequences corresponds 
to a most plausible (i.e., parsimonious) evolutionary history between the genomes 
(see [12, 27]). In the present work we look at the set of simple pseudo-graphs whose 
pressing sequences correspond to the linear extensions of a single poset. Since linear 
extensions can be efficiently sampled asymptotically uniformly, this shows that press- 
ing sequences, and hence the evolutionary histories of the pairs of genomes giving rise 


to said pseudo-graphs, can be sampled near-uniformly. 


Definition 3.2. An ordered simple pseudo-graph, abbreviated OSP-graph, is a simple 
pseudo-graph with a total order on its vertices. In this paper, we will assume that the 
vertices of an OSP-graph are subsets of the positive integers under the usual ordering 
“<” An OSP-graph G is said to be order-pressable if there exists some initial segment 


of V(G) that is a successful pressing sequence. 


Definition 3.3. It was shown in [9] that pressing the vertices of a simple-pseudo- 
graph is essentially equivalent to performing Gaussian elimination with no row swaps 


on its adjacency matrix; therefore, the length of any successful pressing sequence of a 


simple pseudo-graph is the F:-rank of its adjacency matrix. Thus, we define the rank 
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Figure 3.1 Left to right: an OSP-graph G; G1), the result pressing | in 
G; and G?, the result of pressing and then removing vertices 1 and 2 in 
G. Loops are drawn a shaded vertices. 


of a simple pseudo-graph to be the Fo-rank of its adjacency matrix. The rank of a 


simple pseudo-graph on n vertices can vary from 0 (in the case that it is an edgeless 


simple pseudo-graph) to n (such as is the case in Figure 1). We say G is full-rank if 


ni 


its adjacency matrix is invertible over F3. 


Call a matrix M “Cholesky" if there exists an upper-triangular matrix U so that 
M = U'U. In [9] a proof was given that Cholesky decompositions of full-rank, F, 
matrices are unique; in [10] it was shown that for every OSP-graph and adjacency 
matrix A there exists a particular Cholesky decomposition of A that encodes the 


pressing instructions for G. 


Definition 3.4. Let G be OSP-graph with adjacency matrix A (whose rows and 


columns are ordered by the identity permutation). The instructional Cholesky root of 


G (over F 2) is the upper triangular matrix U where for all (i, 7) € [n] x [n], U[, 7] = 1 
if and only if ij € E(Gy_1). In [10] it was shown that U satisfies that U'U = A, 


therefore is a Cholesky decomposition of G. 


The reason this matrix is called “instructional” is that it contains the instructions 
for how vertices affect one another during the corresponding pressing sequence: the 


(i,7) entry is 1 iff pressing 7 flips the state of 7. Since the (instructional) Cholesky 
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matrices are upper-triangular we may also regard U as the adjacency matrix of a 
directed acyclic graph with vertex set {v | v is pressed at some point in the success- 
ful pressing sequence}. Furthermore, the transitive closure of this digraph can be 
considered as a poset. Although it is possible to define these instructional posets for 
less-than-full-rank OSP-graphs, presently we are only concerned with the posets of 
full-rank OSP-graphs. 

We refer to the set of looped vertices in a graph G by L(G) and the set of successful 


pressing sequences for G as }(G). 


Lemma 3.5 ([9], Theorem 9). Let G be a full-rank OSP-graph and o € X(G). Let 


A be the adjacency matrix of G with rows and columns ordered by 0. o € X(G) if 


and only if A has a Cholesky decomposition over F 9. 


Definition 3.6. Let G be a full-rank OSP-graph and o € X(G). Let U be the in- 
structional Cholesky root of A=adj(G), with rows and columns ordered identically 
by o, and D the digraph with vertex set V(G) and adjacency matrix U. The in- 
structional poset of G under a is Poset(G,o) = (V(G), X) where y X x (equivalently 


x = y) if there is an x to y path in D, i.e., Poset(G,c) is the transitive closure of D. 


We say P is generated by G, or equivalently G is a generator of P, if Poset(G,a) = 
P for some o € X(G). If o is the natural order given by G (typically the identity 
permutation) we simply write Poset(G). We denote the set of instructional posets of 


an OSP-graph G by G(G). 


Example 3.7. Let P be a poset on the element set [4] = {1,2,3,4} with cover 


relations 1 > 3, 2 > 3,3 > 4. Then any OSP-graph that generates P must have an 
101 * 


0 11 x 
adjacency matrix A = UTU where U is of the form . It follows that P 


0 Lo 


So 


0 


cS 


0 1 
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Figure 3.2. An order-pressable graph G and the Hasse diagrams 
of the two posets it generates. 


has four generators, as shown below. 
1) (2) 1 3 1 3 1 ae | 3 
b Siw aw wee 
(4) 2 4 2 4 2 4 2 4 
Figure 3.3. The Hasse diagram of P and its four generators. 


We finish this section with two more lemmas from [9] which we will need below. 


Lemma 3.8 ([9], Proposition 1). Let G be an OSP-graph. %(G) 4 0 if and only if 


every component of G containing two or more vertices contains a looped vertex. 


Lemma 3.9 ({9], Theorem 9). Let G be a full-rank OSP-graph and o € 3(G). Let A 
be the adjacency matriz of G with rows and columns ordered by 0. o € X(G) if and 


only if every leading principal minor (over Fy) of A is non-zero. 


3.2. STRUCTURE OF AUTONOMOUS POSETS 


We denote the set of linear extensions of a poset P by LE(P). 
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Lemma 3.10. /f G is a full-rank OSP-graph then LE(P(G,o)) C X(G) for allo € 
Press(G). That 1S, u(G) = Upes(a) LE(P). 


Proof. Let G = (|n],#) be an OSP-graph of rank n ordered by successful pressing 
sequence o. By relabeling G we may assume go is the identity permutation. Let A 
be the adjacency matrix of G (with rows and columns ordered by a) and U be its 
instructional Cholesky root (identically ordered). Let D = ((nl, EB) be the directed 
acyclic graph (aka “DAG”) with adjacency matrix U. Let P = Poset(G) = ([n], Xp) 
and observe that if (a,b) € E then a arent 

Fix a linear extension T = (11,72,...,T) of P. By the previous observation, if 
(aie E then 7; =p 7; and hence 7; must appear before 7; in T = (71, 72,.--, Tn): 


Thus, (7;,7;) € E implies i < 7 € N. By contraposition, we have that 
i > Jj implies (7;,7;) ¢ E. 


Let P be the permutation matrix encoding rt. The previous assertion can be restated 
as 


[PTUP|_=0 for all i> j. 


a,j 


Then V = P’UP is an upper-triangular matrix and 
VTV =(PTUP)T(PTUP) = PTUTUP = PTAP. 


Observe that P? AP is a full-rank symmetric matrix with a Cholesky decomposi- 


tion given by V. It follows from Lemma 3.5 that 7 is a successful pressing sequence 


for G. 


Definition 3.11. We say an OSP-graph G is an autonomous graph if H(G) = 
LE(Poset(G)). We say P is an autonomous poset if there exists an autonomous 
graph G that generates P. That is, if there exists an OSP-graph G such that 
Poset(G, a) = P for some o € &(G) and (G) = LE(P). 
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In our main theorem, we will show that the set of autonomous posets is precisely 
the set of induced N-free and induced bowtie-free posets (referred to in [25] as “V- 


posets”). 


Definition 3.12. For a graph G and a vertex x ¢ V(G) we let x@G be the graph with 
vertex set V(G) U {x}, edge set BGA Cee): and L(x @G) = {x}. Equivalently, 
x ®G is the graph that results from adding a looped vertex x to V(G) and making it 
incident to each looped vertex in G to get an intermediate graph H, then switching 
the state of each edge (including loops and non-loops) in Ny(x) \ {a}. We refer to 
this process as left-appending x to G, we justify this terminology in the following 


observation. 


Observation 3.13. 

Consider OSP-graphs G and H =2@G. Let 7 = (1%,72,..-,;Tr41) € H(A). Since 
L(H) = {x} we have that 7, = x. Furthermore, pressing x switches the state of every 
edge in Ny(x) so H) = G. Thus, the successful pressing sequences of H are exactly 
those resulting from left-appending x to the successful pressing sequences of G. If G 
is order-pressable with instructional Cholesky root U, then x @ G is order-pressable 


and has instructional Cholesky root V that satisfies 
Uli-—1,j—1] ifi,j>2 
Vig] =91 if i =1 andj € L(G) 
0 otherwise. 
Definition 3.14. For a graph G and a vertex x ¢€ V(G) we let G @ xz be the graph 


with vertex set V(G) U {x}, edge set E(G) U {lx | 1 € L(G)}, and 


L(G) if |V(G)| is odd 
L(G@x2x)= 


L(G) U{x} if |V(G)| is even 
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Equivalently, G@ zx is the graph that results from adding a vertex x to V(G), making 
it incident to each looped vertex in G, and, if the resulting graph has an odd number 


of vertices, then we add a loop to x. We refer to this process as right-appending x to 


G. 


Figure 3.4 OSP-graphs z 6 G, G, and G @ a, respectively. 


Recall that the instructional Cholesky root of an OSP-graph is unique. In par- 
ticular, if H is a full-rank graph and V7V is a Cholesky factorization of A = adj(H) 
then V must be the instructional Cholesky root of H; from this we get the following 


observation. 


Observation 3.15. 
If G is order-pressable graph on n vertices and has instructional Cholesky root U 
then G @ x is order-pressable and has instructional Cholesky root V where 
Uli,j] iftjcn 
Vig] = 41 ifj=n+l1 
0 otherwise. 
Lemma 3.16. If G is autonomous then so is x OG. 


Proof. Let H = x@®G. Since £L(H) = {x} we have only one candidate vertex for 


an initial press. Furthermore, by Observation 3.13, H = G. It follows that any 
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pressing sequence must start with 7 and then continue as a pressing sequence for G. 


Therefore, the only instructional poset of H is that of G with a maximum element x 


appended. This demonstrates that H is also autonomous. 
Lemma 3.17. If G is autonomous then so is GQ@ a. 


Proof. If |V(G)| = 1 and G is order-pressable then G is the graph on a single looped 
vertex and G @ x is the graph with one looped vertex, one unlooped vertex and 
an edge between them; both of these graphs are uniquely pressable and therefore 
autonomous. Assume now towards an inductive argument that |V(G)| > 1 and that 
the inductive hypothesis holds for |V(G)| — 1. Let G = ([n],£) and H=G@vz. 
By Observation 3.15, every pressing sequence of G can be extended to a pressing 
sequence for H by appending x to the end of the sequence. Therefore, we need only 
show that |=(H)| = |X(G)| to conclude that H generates only one poset, namely, 
Poset(G) with the addition of a minimal element x. Since Ny(xz) = L(G), the result 
of pressing x (should it be looped) in H would be a loopless graph — by Lemma 3.8 
such a graph cannot be successfully pressed. Thus, every successful pressing sequence 
for H must begin with some element of £(H)\ {x} = L(G). Choose and fix 7 € L(G) 
that is the initial vertex in a successful pressing sequence for H. Assume, by way 
of contradiction, that j is not maximal in Poset(G). It follows that no successful 
pressing sequence for G begins with j, hence (by Lemma 3.8) G) contains a loopless 
component on two or more vertices; call this component C. 


Consider now the result of pressing 7 in H. Since 
Na GQ) OV(C) = L(G) NV(C) = Nala) NV(C), 


we have that every edge from x to V(C) is deleted upon pressing 7 and V(C) is a set 
of unlooped vertices in H. Finally, observe that any vertex that is incident to x in 
H®) must be in a different component than C, as it was in G. It follows that HY) 


contains a non-trivial loopless component, contradicting that 7 was the beginning of a 
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successful pressing sequence. Thus, the initial presses of H are those of G. Observing 


that HY) = GY) © x the result follows from the inductive hypothesis. 


Lemma 3.18. Jf P is an autonomous poset and k is a minimal element, then P —k 


is also an autonomous poset. Furthermore, if 6(G) = {P} then 6(G—k) = {P—k}. 


Proof. Let P is an autonomous poset on n elements. By relabeling, we may assume 
that the elements of P are the integer set [n] = {1,2,...,n}, so that (1,2,...,n) 
is a linear extension of P. By relabeling the minimal elements, we may assume the 
element we remove is n. 

Let G = ((n], ) such that G generates only P. Let A be the adjacency matrix of 
G. By Lemma 3.9 and the fact that G(G) = {P}, for any permutation matrix P we 
have that P? AP has all non-singular leading principal minors (i.e., is LPN) if and 
only if P encodes a linear extension of P. Let A’ denote the (n — 1) x (n—1) leading 


principal submatrix of A. Choose and fix an (n — 1) x (n — 1) permutation matrix 


Pp. 
Suppose P’ A’P’ is LPN. Then 
T 
P' 
a 


is LPN if and only if 


* a 
A! | x 
is invertible. Since A is invertible, we may conclude that if P’A’P’ is 
x |a 
LPN then 
T 
Pg 0) 
A. 
0/1 0} 1 


is LPN. It follows that every successful pressing sequence for a graph G’ with adja- 


cency matrix A’ can be extended to a successful pressing sequence for G by appending 
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n to the end of the sequence. Furthermore, the instructional Cholesky root of A’ is 


the (n — 1) x (n — 1) leading principal submatrix of A; hence G’, the graph whose 


adjacency matrix is A’, generates P — n. 


Lemma 3.19. If P is an autonomous poset and k is a maximal element of P, then 


P —k is also an autonomous poset. 


Proof. Suppose P is autonomous and G is an OSP-graph such that G(G) = {P}. Let 
U be the nxn intructional Cholesky root of G. Then the intructional Cholesky root of 
G is the (n—1) x (n—1) trailing principal submatrix of U. Thus, G™ is a generator 
of P — k. However, every successful pressing sequence of G™ can be left-appended 


by k to obtain a successful pressing sequences for G. Hence, | (e%)| = |X (G)|, so 


that P — k is the only poset generated by G. 


Lemma 3.20. Let P be an autonomous poset on n > 3 elements. If P has a max- 
imum element x and a minimal element z such that x covers z, then any graph G 


that generates only P must satisfy |L(G)| = 1. 


Proof. By assumption that x is maximum we have that P is connected; therefore, if 
y € P \ {x,z}, then x > y and y is incomparable to z. Suppose first that n = 3, 
whence P = ({z,y, z}, X) with x covering both y and z. If G is an OSP-graph that 
generates P then the adjacency matrix A of G must have an instructional Cholesky 
root U encoding the cover relations of P. Hence 
111 111 
U=|0 1 0| andso A=U'U=11 0 1 
00 1 Lo 0 
As the result holds for n = 3, we proceed by induction on n > 4. Choose a 
minimal element y € P \ {x, z}, let P’ = P — y, and let G’ =G-—y. 
By Lemma 3.18, P’ is autonomous and G(G’) = {P’}. Furthermore, P’ has a 


maximum element x and a minimal element z such that x covers z, so we may apply 
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the inductive hypothesis; |£(G’)| = 1, in particular, £(G’) = {x} (since it must be a 
pressable vertex). It follows that L(G) C {x,y}. Assume, by way of contradiction, 
that y € L(G). If zy ¢ E(G) then pressing y would create a looped vertex in 
every component of G™), therefore there is a pressing sequence that begins with y, 
contradicting that P is autonomous. Thus, we must conclude that ry € E(G). Since 
z is a minimal element covered by x, then z is an isolated looped vertex in G) and 
hence Ne(z) = Ne(x). In particular, yz € E(G). 

Let S = Ne(x) \ Ne(y) and T = Ne(y). Assume, towards a contradiction, that 
S #. Observe that sx,sz € E (Go) for all s € S and hence there is a connected 
component in G™ containing x and z (as well as the elements of S$), and z is looped 
in G™. Every other connected component in G) was created by deleting an edge 
between the vertices of 7’ and hence contains an element of T’ which is now looped. 
It follows that G™ can be successfully pressed, which is a contradiction. Thus, we 
may proceed under the assumption that S = 0. 

If v € Na(y) \ Ne(2) then v is looped in G™, ry € E (G), and every other 
connected component in G™) was created by deleting the edge between two unlooped 
vertices and therefore would contain a looped vertex. It follows that Ne(y) = Nea(x), 


therefore x and y can be interchanged in any successful pressing sequence. This 


contradicts that G is autonomous, so we must conclude that y ¢ £(G), as desired. 


Definition 3.21. Let P = (X, x) be a poset. We say (a,b,c, d) is an occurrence of 
the pattern N in P if {a,b,c,d} C X anda>c,a>d, and b> d. We say (a,b,c, d) 
is an induced occurrence of the pattern N in P ifa>c,a>d, b> d and otherwise 
a,b,c and d are pairwise incomparable. 

We say (a,b,c,d) is an occurrence of the pattern bowtie in P if {a,b,c,d} C X 
anda>c,a>d,b>c, and b> d. We say (a,b, c,d) is an induced occurrence of the 
pattern bowtie in P if a > c, a> d, b > d and otherwise a,b,c and d are pairwise 


incomparable. 
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We say P is induced N-free if it contains no induced occurrences of the pattern 
N. Similarly, P is induced bowtie-free if it contains no induced occurrences of the 


pattern bowtie. 


It is worth noting that the literature varies on the definitions of “N-free poset”. In 
our terminology a poset may include an occurrence of the pattern N yet be induced 


N- and bowtie-free. Such an example is the poset P = ([4],{1 > 2>3> 4}). 
Lemma 3.22. Autonomous posets are induced N -free. 


Proof. Let P’ be an autonomous poset and assume towards a contradiction (a, b, y, z) 
is an induced occurrence of the pattern N in P’. Let P = (X, x) be the result of 
iteratively removing maximal and minimal elements from P’ until a,b are the only 
maximal elements and y,z are the only minimal elements. By Lemmas 3.18 and 
3.19, P is an autonomous poset with an induced occurrence of the pattern N, namely 
(a,b,y,z). Observe that if there exists (a’,b’,y’,z') 4 (a,b,y,z) that induces the 
pattern N in P then we may repeat the process of iteratively removing elements 
until only a’,b’, y’, z’ are extremal elements; thus, we proceed under the assumption 
that P has exactly one induced occurrence of the pattern N. 

Choose x € P such that x > y (hence x ¥ y). By assumption that only a 
and b are maximal in P we have that a > x or b> x. Since (a,b, y, z) is an induced 
occurrence of the pattern N we have b ¥ y and hence b ¥ x, therefore a > x. Observe 
that if 7 z then (a,b, x, z) is an induced occurrence of the pattern N, contrary to 
assumption. Thus, x > z (since x # z) and it follows that (x,b,y,z) is an induced 
occurrence of the pattern N implying that 7 = a, therefore a covers y. 

Now choose w € P such that b > w, observe that w #4 a. Since b ¥ y we have 
w xy, hence w > z. Ifa > w > z then (a,b, y, w) is an induced occurrence of the 


pattern N, contrary to assumption. Hence, a > w if and only if w = z. However, if 
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w # z then (a,w,y,z) is an induced occurrence of the pattern N, again contrary to 
assumption. Therefore, w = z and it follows that b covers z. 

By assumption that P is autonomous there exists a graph G that generates only 
P. Fix such a G. Since b € P is maximal, there is a successful pressing sequence 
beginning with b € V(G); thus b € L(G). A sequence o’ = (01,...,0%) is successful 
in G®) exactly when o = (b,01,...,a,) is successful in G. Since G generates an 
autonomous poset then so does G and hence P — b is autonomous. Further P — b 
meets the description of Lemma 3.20 so £ (G) = {a}, therefore L(G) = {a,b} U 
Ng(b). Now observe that if v € Ng(b), then pressing b affects v and hence b > v. It 
follows that Nog(b) = {b, z}, therefore L(G) = {a,b,z}. We proceed to show that z 
can be pressed in G, contradicting that G(G) = {P} 

Suppose first that a € Ng(z). Then Ne(z) \ {b} CL (a) and bv € E (Gc) 
for all v € Ne(z) \ {b}. It follows that any component created by pressing z in G has 
a looped vertex, and hence there is a successful pressing sequence starting with z in 
(G), a contradiction. Thus we must conclude that {a,b,z} C Ne(z). Observe that 
the only elements comparable to y in P are a and y itself. Thus in any successful 
pressing sequence of G, a must be pressed before y and no other vertex affects (or is 
affected by) y. Hence y € £(G) and Ne(y) = Ne(a)\ {y}. Then {a,b, y,z} C Ne(z). 
Since ab, by ¢ E(G) we have that ab, by € E (Gc) and hence a,b and y are path 
connected and y € L (Go), Similarly, if v € Ne(z) \ {a,b, y, z} then bv € E (Go), 
It follows that every non-trivial component created by pressing z in G contains a 


looped vertex, therefore z is the initial press of for some 0 € 4(G), a contradiction. 


Before proceeding, we state the main theorem of [10], which will be used below. 


Theorem 3.23 ([{10], Theorem 1). Let G = (|n], FE) be full rank with instructional 


Cholesky root U. Then G is uniquely pressable (i.e., has exactly one pressing se- 
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quence) if and only if U has columns Cy, ..., C, whose weights (number of nonzero 


entries) are Wi, ..., Wn respectively, satisfying: 


cy = 1, A Neg oe | 


For each j, if Cj = (c1j,Caj,---,Cng)? then 


Gj =0, otherwise 


e Ll=wy, <we < +--+ < Wy. 


w; > 2 implies wi49 > w;, for i € [n — 2]. 


If w; ts odd fori > 1, then w; = for all j > 1. 


For an integer n, let A(n) denote the poset with element set [n] such that n — 2 
covers n and i covers i+ 1 for all i € |[n — 2]. The Hasse diagram of A(n) consist 
of two minimal elements (n — 1 and n) below a chain of length n — 2. Let Gan) 
be the OSP-graph with vertex set V(G) = [n], edge set E(G) = {(i,+1) ie 
[n — 1]} U {(1, 1), (n — 2, n)}. 


Lemma 3.24. A(n) is an autonomous poset and Gain) is the unique graph which 


generates only A(n). 


Proof. Observe that for n = 3 we have only one instructional Cholesky that generates 


A(n); 


aa ames 
= 10 <1 
001 


It follows that the only graph that generates A(3) has adjacency matrix 


‘oe eT 
AS=UAU = 4 6 “4 
hth 


which is the adjacency matrix of G'a(3). 
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For n = 4 we need only consider instructional Cholesky roots of the form: 


1 1 a) 


Ore. abo 
202. cies 20) 
Oe ie 6 ee 


where *1,*2 € {0,1}. A quick check reveals that setting *; = *2 = 0 yields a graph 
with two successful pressing sequences (1,2,3,4) and (1,2,4,3), and otherwise the 
resulting graph has 3 or more successful pressing sequences; hence the claim holds 
form =A, 

We proceed by induction on n > 5. Let G be an OSP-graph that generates 
only A(n). Since A(n) has maximum element 1, we have that 1 € L(G) and GM 
has instructional poset A(n) — 1. But A(n) — 1 is isomorphic to A(n — 1). By the 
inductive hypothesis we have that G™ is isomorphic to G A(a=1) 

Let U be the instructional Cholesky of G under the identity permutation. Let 
A = U'U and let U’ be the (n — 1) x (n — 1) leading principal submatrix of U, 
A’ = U"™U' and G’ = ({n — 1], E’) the graph with adjacency matrix A’. Choose 
and fix 0 € S, such that o(n) = n and let P, be the permutation matrix encoding 
o. Let Pj be the (n — 1) x (n — 1) leading principal submatrix of P, and a’ its 
corresponding permutation. Observe that since G is full-rank then A is invertible. 
Hence, o’ € X(G’) if and only if PZA’P,, is in LPN form, which occurs if and only if 
P? AP, is in LPN form, which in turn occurs if and only if o € X(G). 

Since U(G) = {(1,2,...,n —2,n—1,n), (1,2,...,2 —2,n,n — 1)} we have that 
the only successful pressing sequence of G’ is o’ = (1,2,...,n — 2,n — 1) and hence 
G’ is a uniquely pressable graph (has only one pressing sequence). By Theorem 3.23, 
if U’[1, 2] = 1 then U"[2, 7] = U’|2,4+1] = 1 and hence for 2<i<n—-2 if U[1,7] =1 
then U[2,i] = U[2,i+1] =1. However the intructional Cholesky root of G is the 


(n — 1) x (n — 1) trailing principal minor of U and G™ is isomorphic to Ga(n—1). It 
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follows that U|[2,i+ 1] =0 for all 3 <i <n-—1 thus U[1,i] =0 for all 3 <i<n-1, 
hence [3,n — 1] Ng(1) = 0. Observe that by relabeling n to n — 1 and vice-versa 
we can make the same argument and conclude that n ¢ Ne¢(1), therefore U[1,n] = 0. 


We conclude that G = Gan). 


For an integer n we let X(n) denote the poset with element set [n] so that 1 covers 
3, n — 2 covers n, and 7 covers i +1 for all i € [2,n — 2]. The Hasse diagram of X (n) 
consist of a chain of length n — 4 joining two minimal elements (n — 1 and n) to two 


maximal elements (1 and 2). 
Lemma 3.25. X(n) is not an autonomous poset. 


Proof. Assume, by way of contradiction, that X(n) is an autonomous poset and let 
G be any graph that generates only X(n). Every successful pressing sequence of G 
must begin with 1, 2,3 or 2,1,3. Thus, {1,2} C £(G). Since 3 must be looped after 
pressing 1 and 2, and since the instructional Cholesky root instructs that both 1 and 
2 switch the state of 3 upon being pressed, then 3 € L(G). Observe that X(n)—1 and 
X (n) —2 are isomorphic to A(n—1) and hence G® and G®) are isomorphic to G'a(n—1) 
and hence each have exactly one looped vertex. In particular, £ (a) = 419) for 
{i,j} = {1,2}. Since 1 and 2 are both maximal in X(n) then (1,2) ¢ E(G). It follows 
that Ne(j) = News) for {i, 7} = {1,2}. Therefore, by considering the structure of 
Giain-1), we see Ng(1) \ {1} = Ne(2) \ {2} = {3}; furthermore, L(G) = {1, 2, 3}. 
Consider the result of pressing 3 in G: (1,2), (1,4), (2,4) become edges, 4 becomes 
looped, and every other vertex incident to 3 in G becomes incident to both 1 and 2 
in G®). Thus, there is exactly one component in G®) and it contains a looped vertex 
at 4. By Lemma 3.8 there is a successful pressing sequence in G that begins with 3, 


a contradiction. We conclude that X(n) is not an autonomous poset. 
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Lemma 3.26. Autonomous posets are induced bowtie-free. 


Proof. Let P be an autonomous poset. By Lemma 3.22, P is induced N-free. As- 
sume, towards a contradiction, that (a, b, y, z) is an induced occurrence of the pattern 
bowtie. By iteratively removing maximal and minimal elements, and by application 
of Lemmas 3.18 and 3.19, we may assume a, b, y, and z are the only extremal elements 
of P, and that P does not properly contain another occurrence of the pattern bowtie. 


If the only elements of P are a,b, y, z then 


Sor SS 
Our 2 
So 
Ee Oo 


and hence 


111 £0 
1101 
which has a successful pressing sequence of (4,3, 2,1), contrary to assumption. 
Choose and fix x € P such that x € {a,b,y,z}. Since x is not extremal in P 
we may assume, without loss of generality, thata > x > y. Ifb > x ¥ z then 
(a,b,x,z) induces a bowtie, contrary to assumption. Similarly, if b ~ x > z then 
(x,b,y, z) induces a bowtie. Observe that if b Y « % z then (a,b,x,z) induces an 
N, contradicting Lemma 3.22. Thus we must proceed under the assumption that 
Sa ee 
Observe that the choice of « was arbitrary so any w € P \ {a,b,y, z} must also 
satisfy a > w>yandb>w > z. If x and w are incomparable then (a, b,x, w) and 


(x, w,y, z) induce a smaller bowtie, contrary to assumption. Hence, any two elements 
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in P \ {a,b, y, z} must be comparable, therefore P = X(m) for some m > 5. This 


contradicts Lemma 3.25. 


3.3. MAIN RESULT 


In [21] (and later in [25]) the authors gave a simple description of posets that are both 
induced N-free and induced bowtie-free which we include here as Definition 3.27 and 


Theorem 3.28. 


Definition 3.27. A poset is called a V-poset if it can be generated by beginning with 


the singleton poset and then iteratively applying any of the following three operations: 
(1) a disjoint union, 
(2) adding a new greatest element, 
(3) adding a new least element. 


Theorem 3.28 ({21], Theorem 4.3). A poset is induced N-free and induced bowtie- 


free if and only if it is a V-poset. 


Theorem 3.29. P is autonomous if and only if P is induced N-free and induced 


bowtie-free. 


Proof. By Lemmas 3.22 and 3.26, if P is autonomous then P is induced N-free 
and induced bowtie-free. By Theorem 3.28 it suffices to show that V-posets are 
autonomous. 

A poset on one element is autonomous as it corresponds to the uniquely pressable 
graph on a single looped vertex. We proceed by induction. Let n > 2 and assume 
that all V-posets on n—1 vertices are autonomous. Let P be a V-poset on n vertices. 


If P is the disjoint union of multiple posets then each of its connected subposets is a 
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smaller V-poset. By inductive hypothesis for each connected subposet there is a graph 
that generates it and has only the pressing sequences dictated by said subposet. It 
follows that in this case P is autonomous as well. Suppose now that P is connected. 
It then follows that P has a unique maximal or a unique minimal element. Let P — x 
be the result of removing a unique maximal or minimal element from P. Observe 
that P — x is a V-poset and thus by induction is autonomous; let H be a graph such 


that G(H) = {P — x}. By Lemmas 3.13 and 3.15, 7 @ H or H @z generates only P 


and therefore is autonomous. 


3.4 V-POSET RECOGNITION 


For a poset P we let np and ep denote the number of vertices and edges in the 
Hasse diagram of the poset, respectively. We let hp denote the sum of the heights of 
components of P (the height of a poset is the length of its longest chain), cp denote 
the number of components of P, and Mp and mp denote the number of maximal and 


minimal elements in P, respectively. 


Lemma 3.30. If P is a V-poset then 
ep = 2np+cp — Mp — mp — hp < 2np — 2 


Proof. We show that ep = 2np+1—Mp—mp-—hp for a connected poset; the equality 
above follows by summing over components, and the inequality is immediate. Observe 
that if np = 1 then P is a poset one element and hence (2np+1)—(Mp+mp+hp) = 
0 = ep. Assume towards an inductive argument that np > 2. Since P is connected 
it must have a unique minimal or maximal element, say x, which we assume will be 
maximal (as the argument is identical for a minimal element). Let Q = P — {zx}. 


Then, by applying the inductive hypothesis to Q, 


ep — Mg = €g = 2ng+1-—Mg-—mg-—hog 
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ep = 2n9 +1—Mg — ho = 2(np —1) +1—mp — (hp — 1) 


ep = 2np — mp —hp 


By noting that Mp = 1, we have our result. 


We now give a different edge count that uses width (referred to as wp in the 
statement) instead of heights. While both of these edge counts are necessary for the 


property of being a Y-poset, even when taken together, they are not sufficient. 
Lemma 3.31. If P is a V-poset then 
ep = Np +wp—- Mp —mp 


Proof. As in the previous proof, we show that ep = np + wp — Mp — mp for a 
connected poset; the equality above follows by summing over components since the 
width of a disconnected poset is the sum of the width of its connected components 
(i.e. the length of a maximal antichain). Observe that if np = 1 then P is a poset one 
element and np + wp — Mp — mp = 0 = ep. Assume towards an inductive argument 
that np > 2. Since P is connected it must have a unique minimal or maximal element, 
say x, which we assume will be maximal (as the argument is identical for a minimal 


element). Let Q = P — {x}. Then, by applying the inductive hypothesis to Q, 


ep = eg + Ma = (ng + wa — Mg-— ma) + Mo 


=NngtWog-Mg=np—l+wp-—mp=np-— Mpt+wp- mp. 


We propose an algorithm for the recognition of autonomous posets that operates 
on an arbitrary directed acyclic graph whose transitive closure is the poset in question. 
As a subroutine, we employ an algorithm found in [29] that detects if a directed 
acyclic graph contains an induced copy of the pattern N and, if the input is found 


to be induced N-free, it also returns the transitive reduction of the input. The 
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aforementioned subroutine is guaranteed to run in O(|V| + |£]|). Observe that by the 
proof of Lemma 3.26, in order to determine if an induced N-free poset is a V-poset 
we need only to verify that its transitive-reduction does not contain a sub-DAG that 
is isomorphic to ((4], {(1, 3), (2,3), (1, 4), (2,4)}) (as done in Subroutine 2) and does 
not contain sub-DAG whose transitive closure (interpreted as a poset) is isomorphic 
to X(n), (n > 5). 

Lemma 3.30 shows that if we present the poset by the transitively-reduced directed 
acyclic graph with cover relations as edges then the run-time is O(|V|). Observe that 
in Subroutines 2 and 3 each edge is traversed at most twice, hence these algorithms 
have run-time O(|V|+ ||). Thus the presented algorithm has the same run-time as 


Subroutine 1. 


Algorithm 1 
1: input: a directed acyclic graph D. 
2: output: true or false. True if the transitive closure of D is a V-poset, False 
otherwise. 
3: Bool + true 
4: if IsSeriesParallel(D)|[Bool]=False then 
5: Bool + false 
6: else 
7: D+ IsSeriesParallel(D)|DAG} 


8: if IsBowtieFree(D)= false then 


9: Bool + false 

10: else 

iM if ClosurelsVPoset(D) = false then 
12: Bool + false 


13: return Bool 
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Subroutine 1: IsSeriesParallel() 
1: input: a directed acyclic graph D. 
2: output: (Bool, DAG). Bool= true when D has a series-parallel decomposition 
and Bool= false otherwise, and DAG is the transitive reduction of D. 
3: Algorithm found in [29] 
4: return (Bool, DAG) 


Subroutine 2: IsBowtieFree() 


1: input: an induced N-free, transitively reduced directed acyclic graph D. 

2: output: true or false. False if some induced subgraph of D is isomorphic to the 
bowtie digraph ({a, b,c, d}, {(a,c), (a, d), (b,c), (b, d)}), True otherwise. 

3: Bool < true, Current «+ @, Parents «+ Q, Visited «+ @ 

4: for v € V(D) do 

5: if OutDegree(v) = 0 then 

6: Current.Add(v) 

7: while Current4 @ do 

8: for v € Current do 

9: for u € InNeighborhood(v) do 

10: Parents. Add(w) 


11: for v € Parents do 


12: if OutDegree(v) > 1 then 

13: for u € OutNeighborhood(v) do 

14: if InDegree(u) > 1 then 

1b: Bool + false (Break while loop) 
16: Visited. Add(w) 


17: for v € Current do 


18: Visited. Add(v) 
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19: 


20: 


21: 


22: 


23: 


24: 


25: 


Current + 0) 
for v €Parents do 
if vu ¢ Visited then 
Current. Add(v) 
Visited. Add(v) 
Parents + @ 


return Bool 


Subroutine 3: ClosureIsVPoset() 


1: 


input: an induced N-free, induced bowtie-free, transitively reduced directed 


acyclic graph D. 


: output: true or false. True if the transitive closure of D is a V-poset, False 


otherwise. 


: Bool + true, Current + 0, Parents «+ 0, Visited + @, Multiple + @ 


: for v € V(D) do 


if OutDegree(v) = 0 then 
Current. Add(v) 


: while Current 4 @ do 


for v € Current do 
if v € Multiple and InDegree(v) > 1 then 
Bool + false (Break while loop) 
for u € InNeighborhood(v) do 
Parents. Add(w) 
if v in Multiple then 
Multiple. Add(w) 
for v € Parents do 
if OutDegree(v) > 1 then 
Multiple.Add(v) 
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18: 


19: 


20: 


21: 


22: 


23: 


24: 


for v € Current do 
Visited. Add(v) 
Current + 0) 
for v €Parents do 
if vu ¢ Visited then 
Current. Add(v) 


Parents + @ 


25: return Bool 
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CHAPTER 4 
CHOLESKY ROOTS 
4.1 INTRODUCTION 


Over the complex field C, a square matrix M is said to have a Cholesky factorization 


if there exists an upper-triangular matrix U so that U*U = M, where * denotes 


the conjugate transpose of a matrix (or simply transpose when restricting to R). 


For a prime power qg, we say M has a Cholesky factorization if there exists an upper- 
triangular matrix U so that U'U = M. Observe that not all matrices have a Cholesky 
factorization. For example, only symmetric matrices have Cholesky factorizations 
since 
hy ai SO 

For a matrix with complex entries there are a multitude of equivalent characteriza- 
tions that determine if a matrix allows a Cholesky factorization (see, e.g. [22]). One 
particular example of this is the notion of positive-definiteness which, with some care, 
can be extended to some finite fields. For a wonderful survey, and some surprising 
results, on positive-definiteness over finite fields see [19]. In [9], the authors describe 
a surprising connection between successful pressing sequences of an ordered simple 
pseudo-graph (OSP-graph) and Cholesky factorizations of its ordered adjacency ma- 
trices. In particular, they argue that an OSP-graph G with vertex set {v1, v2,...,Un} 
can be pressed in order vj, v9,...,Un exactly when the adjacency matrix A of G, 
with rows and columns ordered by v1, v2,...,Un, has a Cholesky factorization. This 


equivalency only holds for OSP graphs satisfying that each vertex appears in (each) 
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successful pressing sequence. In Chapter 2, we define a special instructional Cholesky 


factorizations that extends the above equivalency to OSP graphs whose pressing se- 


quences do not include every vertex. 


In the following sections we explore how many distinct Cholesky factorizations 


exist for a matrix with entries from Fy. 


4.2 A BIJECTIVE ARGUMENT 


For all positive integer n, we let 1, and 0, denote the n x n, 


F, multiplicative and 


additive identity matrices (respectively). For r < n, we let U,(r) be the set of n x n, 


rank r, upper-triangular matrices with entries from F2. For n > 1 and r < n we 


define 
A,r) S40 et). ||US 14 
Yalr) = {U €U,(r)|U? = 0} 
Zilr) = {U eU,(r) |UTU =0,} 
and 


Xn= DL lAn(r) Ya= Do alr), 42n= dL [Fnl(r)| 


O<r<n O0<r<n O0<r<n 


Observation 4.1. For all n > 1, Xp, = Yn. 


Proof. Let y : F3*" — F5*” by y(A) = A+1,. Then 


Theorem 4.2. For alln >1 andO0<r<n: 


[Yn(r)] = |Fn(r)| 
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Proof. Observe that 

v1(0) = {0} } = 2:00), 
and Vi(r) = Z1(r) = @ for all r 4 0. We proceed by induction. Let n > 1 and assume 
that |V,»-1(7)| = |Zn-1(r)| for all r <n — 1. Choose and fix a rank r, n x n upper- 
triangular matrix B. Observe that by Sylvester’s rank inequality Y,,(n) = Z,(n) = 0, 
so we may proceed with the assumption that r < n. Let B’ be the n-—1xn-1 


principal submatrix of B. 


Vv B? | Blv +bv 


07 |b or | b or | b? 
Then B € ),, if and only if b = 0 and B’v = O and B’ € J,_1. However B’v = 0 
if and only if v € Null(B’), the null space of B’. If B’ € Y,_1 then the column 
space of B’, Col(B’), must be a subset of Null(B’). It follows that if B € Y, then 
v € Col(B’) or v € Null(B’) \ Col(B’). 


It follows that for each r: 
IVa (r)| = |Yw-1(r)| - 2” + |[Yn-i(r — 1)] - (gene) = a) 


Yn (r)| = [Yw—1(r)| 2" + |Yw-a(r — 1) - (Ce = aa 


Choose and fix a rank r, N x N upper-triangular matrix C. Let C’ be the N—1x N-1 


principal submatrix of C. 


CT (8) C' Ww CTC’ C Tw 


CCc= 


Cc Oc wiC' |wiwt 


C € Zy if and only if w7.w+ c? = 0 and w’C’ = O and C’ € Zy_,. Equivalently 
C € Zy if and only if C’ € Zy_, and 
CIT 


0 w 


Cw= 


1 | Cc 
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This occurs if and only if w € Null(C). Observe that if c = 0 then w € Null(C) 


exactly when w € Row(C’)N Null(C) or w € Null(C) \ Row(C’). On the other hand 


if c= 1 then w € Null(C) exactly when w € Null(C) \ Row(C’) . 


It follows that for each r: 
IZn(r)| = |Zn_1(r)| - 2” + |Zn_1(r — 1)]- (oaaCruie)) 7 a) 


|Znx(r)| = [Znv-a(r)| 2” + [Zn-alr = DI (24 = 24) 


Theorem 4.3. 
Y= “ es) - aire) | gr ese —F, if n is even 
: (aes a iu, asa Di Ia DSI 23 if n is odd 


Proof. Theorem 1 in [13]. 


Corollary 4.4. 


ae . wat sy) — (nyr sy) | 2, ifn is even 


n+l n n 2/44(n 5294 : 3 
* (eer, = re gin" /4+(n—-1)/2-37" 25, ifn is odd 
Corollary 4.5. For alln > 1, Yn(r) = Z,(r) = 0 whenever r > n/2. 


Proof. If A € Y,z(r) then Col(U) Cc Null(U) where the inclusion is strict since 


[0,...,0,1]7 € Null(U)\Col(U). That is r <n —r. Since there is a rank-preserving 


bijection between Y,,(r) and Z,,(r) the result holds. 


Corollary 4.6. Let A be a symmetric, rank r matrix with entries from Fz. For each 
k € |n], let Ay be the leading principal submatrix of A. If A satistfies that det(A;) = 1 


if and only ifk <r (a.k.a. A is in leading principal minors non-negative form) then 
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the number of distinct Cholesky factorizations for A is 


“ Gee) = Cosa) QA es, if n is even 


» [ae = acti acn)| g(n-1)?/4+(n-1)/2-3?-25 i mis odd 


where n is the corank of A. 
Proof. Let r be the rank of A and let A; be the principal r x r submatrix of A: 
Ai 


Aj2 Aj, | A1,2 


Ag» 


T 
Aoi Ao» Aj» 


If B?B = A is a Cholesky factorization of A then 


T T 
Bi 2 By Pia By Pie 


Bry 0 Bia 


BTB= 


Bt) Be sd| |. 0: || Bye || © |sBEBig| BEBE BT Bes 


But A has an instructional Cholesky of the form 


T T 
Vie Via | ViaMie 


Vii 0 Via 


ViIV= 


vz,;o}] 0 | o VEGA VAs 


By uniqueness of Cholesky decompositions of full-rank matrices over GF'(2) we have 


that By; = V1. Then by invertibility we have 


Bia = (Vi) BE Bia = (VR) ViiVia = Viz 


and hence 
ViQMi2 = BL Biz + By Bo2 > By» Bo2 =0 
Then 
eile = Vi, 0 Via Via = A 
eternal eater 


if and only if CZ C22 =0 


ws 


4.3. ASYMPTOTIC BEHAVIOR 


2n 
Here we will look at the asymptotic behavior of Z. Let C, = (a) denote the n'* 
Catalan number, let (x), = x(a — 1)--- (a —k +1) denote the descending factorial, 


and let 


Fis _ - 2n 67 a 1 Q-5(35+1) = > 2n (2 + ~) Q-37°-5 
7 NE 37) nae pel ae A n+1 


Then 
n+1)3c44 posses 
Bey | SS ele, J 1) (67 +1)2°9 3 
Bon Oe ieee erat) 
and 
(n+1)3j41 ntl ate ae 
(n+ 374+ 1)aj41. nt+3j7+1 n+1 n+1 
However 
ee = 4 ~3j ae 3j +1 (-1)'( 37 i 
‘fal n+1 =5 t n+1 
s ne (1 (92) a1 BOD OF 
ram n+l n+l 7 Cae al 
Hence 
7.671 
|Fon| < S> 997(67 + 1)27 Sy 
nel ui<m/3 cae 
Thus, 
ae _ SS 2n _ 2n gn? —37?—j 
j n— 37 woop as 
Zon 


= 3 2n 67 +1 9-37?—j 
Qn? j n—37/ \n+374+1 


Lon 67 +1 372-4 
= ee eee 
gn? we sacs 
Zon 2. Eon 
me =, S- (6) + 1) g-8r 4 + ae) 
j n 
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Using Wolfram Alpha yields 


and 


It follows that 


Moreover, 


lim > (6j+1)2°°7-9 = 0.17755 


n> 
lj|<n/3 
. Lon 7.671 
lim —— < lim _ 


Zon 0.17755 - 2" C., 


MoV 7 os ONT 58 0 Ce 
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CHAPTER 5 


UNIFORM SAMPLING ALGORITHM 


5.1 INTRODUCTION 


In 1969, Jesse MacWilliams published a paper counting the order of the subgroup 
of matrices U in GL(n,q) that satisfy UU? = I. In order to do so, he first gave 
a recursive algorithm that counts the number of symmetric matrices with entries in 
GF‘(q) of size t x t and rank r. By reversing this algorithm we are able to construct 
a method for uniformly and efficiently sampling symmetric matrices of a given size 
and rank over a finite field. In the context of simple pseudo graphs his means that we 
can randomly construct a simple pseudo-graph on t vertices whose successful pressing 


sequences will be of length r. 


5.2 ALGORITHM FOR SAMPLING 


We begin with some notation that will be useful in understanding the construction 
of sampling algorithm. 


Notation: 
q denotes a prime power, 
GF‘(q) is the finite field of g elements, 
GL(n,q) the group of n x n invertible matrices with entries in GFj, 


O(n, q) the subgroup of matrices U in GL(n,q) that satisfy UU? = J, 
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S(t,r,q) the group of symmetric matrices with entries in GF'(q) of size t x t 


and rank r, 


N(t,r,q) the number of symmetric matrices with entries in GF'(q) of size t x t 


and rank r. 


The following key lemma was proven in [24]. Nevertheless, we present a brief 


proof here since it is instrumental in constructing the sampling algorithm. 


Lemma 5.1. Let A be a symmetric t x t matrix of rank r with entries in GF(q) and 
let N;(A) denote the number of symmetric (t+ 1) x (t+1) matrices of rank r+i with 


entries in GF'(q) that contain A as a leading principal submatriz. 


ie _ qt, if i —9 
(q—1)q’, ifi=1 
N;(A) = 

q', ifi=0 

0, otherwise. 
Aly? 

Proof. Fix A and consider a matrix of the form where y is an 1 x t row 

Y | Yo 


vector with entries in GL(q) and yo € GL(q). If y is not linearly dependent on the 
rows of A then the resulting matrix has rank r + 2. There are g' — q” such choices 
for y. Since we are free to choose any entry for yo this yields No(A) = q(q' — q’). 
Suppose now that y is a linear dependent on A, that is y = nA for some 1 x t vector 
n. If yo = ny’ then the resulting matrix has rank r, otherwise if yy 4 ny? then the 
resulting matrix has rank r +1. Since there are q’ choices for y in span(A) we have 


N,(A) = (¢ — 1)q” and No(A) = @’. 


Observe that from this we obtain 
N(t+1,7r,q) = (gt! -@") N(t,r — 2,4) + @—- Da N(t,r -1gt+¢ Nr) 
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Fix t,g,r where t > 0, g a prime power, and 0 < r < t. Assume that for 
all r’ < r we can sample S(t,r’,q) uniformly at random. Consider the following 


sampling algorithm for S(t+1,r,q). Let 


Tr N(t,r,q) j— 
GT Nt”? i=0 


— r—1 N(t,r-1, ) a 
Pi (q—1)q¢ NE CLRA i=1 


tay ght rhe i-2 


(q N(t+1,r,q)’ 


and observe that 


i _INGr@ +(a- Vg NG -La+ gg NG = 2.9 _ | 
L N(t+ ie, q) , 


Select 7 with probability p; then take a sample from S(t,r—i,q) uniformly at random, 
call this matrix A. If 7 = 2, choose y uniformly at random from the q‘ — q” vectors 
that are linearly independent of the rows of A, choose yo uniformly at random from 
{0,1,...,q¢—1}. Otherwise, if i 4 2 choose, uniformly at random, a vector n with 
entries in GF'(q) and set y = nA. Furthermore, if i = 0 set yo = ny? and if i = 1 


choose yo uniformly at random from {0,1,...,q— 1} \ {ny7}. 
P 
Let A’ = , by the proof of Lemma 5.1 A’ is a uniform selection among 
Y | Yo 
the matrices in S(t + 1,r,q) who contain A as a leading principal submatrix. 


Fix positive integer n and integer 0 < r <n. If r = 0 there is only one matrix to 
consider. Otherwise, choose i € {0, 1,2} with probability p; := p;(t,r,q) (above). Set 


k = 1. While n > 1, repeat the following process: 
e set 2, =1 
e replace t with t — 1 
e replace r with r —7 


e select and replace i € {0, 1,2} with probability p; 
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e replace k with k+1 


Clearly this process ends since at each step t is reduced. If at the end of the process 
r = 0, set Ap = [0]. Otherwise, choose a uniformly at random from {1,2,...,q— 1} 
and set Ap = [a]. for each i € k, construct Aj,; from A; by uniformly increasing the 
rank of A; by 7,_;. This results in a t x t matrix A,_, with rank r. 

Choose and fix r > 0 and t > 0. Select A uniformly at random from S(t+1,r,q). 
Observe that if t = 0 the algorithm selects with A with probability a Assume t > 1 
and that the algorithm selects uniformly for all values up to t — 1. The probability 


that the algorithm selects A is 
P(A) = P(A|i = 0)P(@ = 0) + P(Alé = 1)P( = 1) + P(Ali = 2)PC = 2) 


Let B be the t x t leading principal submatrix of A. 


1 1 


P(A|i) = P(B € S(t,r —1,q)) 


where N;(B) is the number of symmetric t+ 1 x t+ 1 matrices of rank (r — 7) +i 


that contain B as a leading principal submatrix. 


NG e159) 


P(B —14,q)) = 
(B € S(t,r —i,q)) N(t,r—2,q) + N(t,r —1,q) + N(t,7,9) 


so 
a N(t,r — i,q) i i 
P(A = é J ni 
( ) » NGr=t NEF ene Grae 
_ i Nit,r,qpo , N(tjr—-lan , N(t,r—2,¢)po 
X \N(t,r,q)No(B) © N(tr—1,q)Ni(B) " NG — 2,4)No(B) 
= 1 Nit,7,q) | ING = 1@) \ N(t,r — 2,q) _ it 
~ X\NG+17q) °° NE+1L7q ° N(t+1rq)}]  N(t+1,7,9) 


where X = N(t,r — 2,q)+ N(t,r-—1,¢) + N(t7,¢) 


Below, we present the algorithm in more detail. 
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Algorithm 1: NewRandomMatrix 
1: input: An integer tuple: (n,r,¢) 
2: output: An x n, rank r, symmetric matrix with entries in F, 
3: if n =1 then 


4: if r=O0 then 


Bi M= 0] 

6: else 

7: x «-RandomNonZero(F,) 
8: M= c 

9: else 


10: for 7 in {0,1,2} do 

11: pi < ProbFrom(n, r, 1) 

12: Choose an element of {0,1,2} with probability (p1, po, p3) 
13: M «— NewRandomMatrix(n — 1,r — j,q) 

14:  M + RanklIncrease(M,j) 


15: return M 


Algorithm 2: RankIncrease 


1: input: An n x n, rank r, symmetric matrix M with entries from F, and an 


integer 7 € {0,1,2} 


2: output: An (n+ 1) x (n+1), rank r+ J, symmetric matrix with entries from F, 
3: if j=2 then 
4: V < NotInSpanVector(M) 


4 


5: u« RandomElement(F,) 
v {Vt 


6: return 


V| M 
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ts 


8: 


9: 


10: 


11: 


12: 


13: 


14: 


15: 


else 
V < InSpanVector(V/) 
i yaa r Vad 
if j=0 then 
for i € [n]: r; <— RandomElement(F,) 
Ue" 


else 


Choose v uniformily at random from F, \ {x} 


v|Vr 


return 


V|M 


Algorithm 3: NotInSpanVector 


: input: An n x n matrix M of rank k with entries in F, 


: output: Ann x 1 vector in the complement of span(V/) 

: L’ + FindOutBasis(/) 

: Set = n-—k and choose a random integer r from [1, q‘ — 1] 
: Cast r to a q-ary vector R 

. R< FIELD(R) 

: V + InSpanVector(V/) 


: return [/-R+V 


Algorithm 4: InSpanVector 


1: 


2: 


3: 


input: An n x n matrix M of rank k with entries in F, 
output: Ann x 1 vector in the span of MW 


L + FindInBasis(M) 
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4: for 7 € [1, k] : 7; — RandomElement(F,) 
iy 
5: V= T, To ++: Tk 


6: return L-V 


Algorithm 5: FindOutBasis 

1: input: a matrix MV 

2: output: a basis for the complement of 
3: B <FindInBasis(/) 

4: find C to complete B to an n x n basis 


5: return C 


Algorithm 6: FindInBasis 
1: input: a matrix MV 
2: output: a basis for M 
3: use in place Gaussian elimination to find column basis to find basis B 


4: return B 


Algorithm 7: ProbFrom(n,r,i) 
1: input: Integer tuple (n,r, 7) 
2: output: A probability p 
3: Denom + Numb(n — 1,7r,0)+Numb(n — 1,r — 1,1)+ Numb(n — 1,r — 2, 2) 


4: return Numb(n — 1,r —i,i)/Denom 
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Algorithm 8: Numb(n,r, 7) 


1: 


2: 


input: Integer tuple (n,r,7) 


output: The number of symmetric (n+ 1) x (n+ 1) matrices of rank r +i with 


entries in F, that can be generated by appending the same vector as a row and 


column, along with a corner entry, to any symmetric n x n rank r matrix with 


entries in F, 


: if i = 0 then 


return NumbSymm(n,r,q) q” 


: else if i= 1 then 


return NumbSymm(n,r,q) (q¢ — 1)q" 


: else if i = 2 then 


return NumbSymm(n,r,q) (q”*1 — q"*!) 


Algorithm 9: NumbSymm 


1: 


input: An integer tuple (n,r,q) where n > r > 0 and q is a prime power 


: output: the number of symmetric rank rn x n matrices over 
ahs PA 
faa? =) 


CG Nee a | 


: return hes f(i)) near g(2)) 


Algorithm 10: RandomNonZero 


1: 


2: 


3: 


input: q 


output: a random element of F, \ {0} 


return Random element from F, \ {0} 
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CHAPTER 6 


THE PRESSING SPACE OF A GRAPH 


6.1 INTRODUCTION 


In 1999, Hannenhali and Pevzner demonstrated that a shortest path of reversal ed- 
its from one permutation to another can be determined in polynomial-time. The 
polynomial-time argument is justified by observing that a set of permutations can 
be converted into a simple pseudo-graph in polynomial-time (see the introduction of 
this thesis) and then the length of such a conversion path is simply the rank (over 
F,) of the adjacency matrix of the graph (which can be done in sub-cubic time). On 
the other hand, this does not answer how many shortest path of reversal edits exist. 


Thus, we propose the following question: 


Question 6.1. What is the computational complexity of determining the number of 


successful pressing sequences of a simple pseudo-graph? 


For a full-rank, ordered simple psuedo-graph G on n vertices, there are n! per- 
mutations to consider. Attempting to press G in order o is equivalent to performing 
Gaussian elimination without pivoting on the adjacency matrix of G, this can be done 
in polynomial-time. Since the average graph has oa successful pressing sequences, the 


exhaustive method is not a practical computation technique for determining |G(G)|. 


6.2. BLOCK PRESSING AND EXPONENTIAL COMPLEXITY 


In what follows we will give a method that demonstrate that |G(G)| can be determined 


in exponential-time. For simplicity of the argument we (re)introduce some notation 
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here. 

A pressable graph is a simple pseudo-graph with a looped vertex in every non- 
trivial component. A vertex in a pressable graph is said to be stuck if it is looped 
but pressing the vertex creates a loopless graph with at least one edge. A vertex in 
a pressable graph is said to be eventually stuck if it is looped but pressing the vertex 
creates a loopless component with at least one edge. A looped vertex that is not 
eventually stuck is referred to as a pressable vertex. In {9], the authors showed that 
if G is pressable graph then we can always find a pressable vertex by creating a set 
X of looped vertices with the fewest possible number of neighboring looped vertices 
and then selecting « € X such that deg(x) is maximal in X. 

Let V’ C V(G) and let G’ = G[V’] the induced subgraph of G. If G’ is a pressable 
graph then we can find a pressing sequence by iteratively finding a pressable vertex 
and pressing it. We say V’ is block pressable if there exists a successful pressing 
sequence o’ of G’ that can be extended into a successful pressing sequence in G. 
That is, if pressing V’ in the order given by o’ does not create any (eventually) stuck 
vertices in G. For the following lemma we let G,, denote the result of pressing the 


vertices in appearing in o’ in G. 


Lemma 6.2. Let G = (V, E) be a simple pseudo-graph and let V' C V and suppose 
G’ is block pressable. Then G, = G,-for all o',7’ € G(G’). 


Proof. Let G = (V, FE) be an ordered simple pseudo-graph and let G’ = (V’, E’) be 
a block pressable subgraph of G. By re-indexing, we may assume that both G and 
G’ are identity pressable. Let A be the adjacency matrix of G and A’ the adjacency 
matrix of G’. Recall that a gaph is order pressable if and only if its ordered adjacency 
matrix is in leading principal minors non-negative form ({9]). 


Then 


and both A and A’ are in leading principal minor form. Let 7’ be any successful 


pressing sequence of G’ and let P,, be the permutation matrix encoding r’. Then 
PT AP, 


is in leading principal minor form. Consider 


. feo] fa 
PPAP = 
Pole] [a 


B|| Pp, o| | ptare, Pre | 
el[o lr] farm [e | 


Then the columns of P4B are linearly dependent on P4A’P,, exactly when the 


columns of B are linearly dependent on A’. It follows that r’ can be extended to 


a successful pressing sequence of G. 


Example 6.3. The vertices in G’ can be be pressed in two orders: (v1,v3, v2) and 


(v3, U1, V2). The resulting graph is the same. 


Figure 6.1 Block Pressing 


Theorem 6.4. Let G be a full-rank, ordered simple pseudo-graph. Then |G(G)| can 


be determined in exponential-time. 


Proof. Let G = (|n], £) and let B,, be the boolean lattice on the set {1,2,...,n}. Let 


H = (V(H), E(H)) be the Hasse diagram of B,,. 
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We define two functions iteratively, one on the elements of H and one on its cover 
relations. f : E(H) > Z and g: V(H) \ @ > Z defined iteratively by 
1, if x is pressable in G 


F({O, {2} }) = and g({x}) = FLO, {2}}) 


0, otherwise 


for each x € [n]. 
For each S$ C [n] and y ¢ S let 
1, if g(S) £0 and y is pressable in G after pressing S$ 
FS, SU {y}}) = 


0, otherwise 


and for each T C |[n] 


gL) = DFAT \ {t}, TH) 9 \ (th) 


teT 


Observe that when g(T) 4 0 it records the number of successful way to press the 
vertices of the induced graph G/T]. Hence g({n]) = |6(G)|. The number of points 
considered by f is n2"-' and the number of points considered by g is 2”. Each 


function evaluation consisted of addition (linear time) and block pressing which is 


polynomial-time. Therefore, the time complexity of computing g([n]) is O(2"). 


6.3. THE PRESSING GAME CONJECTURE 


A sequence 7 = 7,...,7% is a subsequence of permutation o = (01,...,0n) if 
(Ci Cassia) = (Mig. 7e) for some lS <i de. <ee ee Sm. These 
quence 7 is a common subsequence to o and T if it is a subsequence of both o and T. 


We let Ics(o,7) denote the maximum length of such a subsequence, that is 
les(a;7) = max ((71 : 7 is a common subsequence to o and T) 


Given two permutations, o and T, the edit distance between o and 7, denoted d(a,7T), 


is the minimum number of entries that must be removed and replaced (in any position 
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Figure 6.2 Computing the pressing number in exponential time 


and order) from o to achieve 7. A simpler formulation is 
d(a,T) =n — lcs(a,T) 


where |o| = n = |r|. 

Let G = (|n], £) be an ordered simple pseudo-graph. The pressing space of G, 
denoted G(G), is the set of all successful pressing sequences of G. From this we 
create a family of metagraphs {M;,(G)} where V(M;,(G)) = G(G) and E(M;(G)) = 
{{o,7} | d(o,7) < k}. 


In [4] we find the following conjecture. 


Conjecture 6.5. Every successful pressing sequence can be reached from every other 


one by a sequence of edits that involve at most four deletions or insertions. 


Using the terminology establish above the conjecture can be rephrased as M4(G) 
is a connected graph for all ordered simple pseudo-graphs G. The conjecture has been 
shown to be true on very specific types of graphs (see [4]), but remains unresolved in 


general. The following example demonstrates that M3(G) is not connected. 
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Example 6.6. Let G = (V,E) be a simple pseudo graph with V = {1,2,3,4,5,6} 
and F = {(1, 1), (1, 2), (1, 5)(2, 3), (3, 4), (3, 6), (5, 5), (5,6)} where we use parenthesis 
to avoid ambiguity on the looped vertices. With some effort we determine that 


6(G) = {(1, 2,3, 4, 5, 6), (5, 6, 3,4, 1, 2)}. However, 
a((l. 2,34. 5,6))(5,6,3,4,1,2)) =4 
since any longest increasing subsequence in (5,6, 3,4, 1,2) is of length 2. 


This previous example demonstrates that if the pressing game conjecture is true 
(at distance 4) then it is best possible. While the relatively small size of the previous 
example may cast some doubt on the conjecture, it has held up (so far) to all tests. 
In the next section, we (joint work with Joshua Cooper and Peter Gartland) present 


a weaker conjecture and a proposed method to resolve it. 


6.4 ‘THE WEAK PRESSING GAME CONJECTURE 


Conjecture 6.7 (WPGC). Every successful pressing sequence can be reached from 
another successful pressing sequence by a sequence of edits that involve at most four 


deletions or insertions. 


The previous conjecture is implied by the pressing game conjecture and can be 
restate as such: 


WPGC: M,(G) has no isolated points 


In what follows, we give an argument as to why Conjecture 6.7 should be true. 

In Chapter 2, we describe the set of uniquely pressable graphs by characterizing 
their Cholesky roots. Below is the characterization again, but using some more 
colloquial language. 

An nxn upper-triangular matrix U is the Cholesky matrix of a uniquely pressable 


graph provided 
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i. U has ones along the diagonal and super-diagonal; 
ii. in the upper-triangular portion of a matrix, a one is never above a zero; 


iii. if column j contains J > 2 ones, then column j +1 contains at least J ones (for 


all 7 < n—1), and column j + 2 contains at least J + 2 ones, (if 7 <n — 2); 


iv. if column 7 > 1 has an odd number of ones, then column k has only ones above 


the diagonal for each 7 <k <n. 


Definition 6.8. We say that an n x n upper-triangular matrix U is minimally non- 
uniquely pressable (hereafter referred to as MNUP(n)) if it is not the Cholesky matrix 
of a uniquely pressable graph, but every one of its principal submatrices is. We say 
G is a minimally non-uniquely pressable graph, or simple MNUP graph, if its identity 


ordered instructional Cholesky is minimally non-uniquely pressable. 


One can easily verify that there are no MNUP graphs on 3 vertices. Let M be 
a MNUP(n) matrix for some n > 3. Let U and V be the (n — 1) x (n — 1) leading 
and trailing principal submatrices of M, respectively. Since U is uniquely pressable 
we have that the first n — 1 columns of M satisfy uniquely pressable laws 2, 72, 277, av. 
For simplicity of argument, for any matrix A let a(A) denote the entry in the first 
row and last column of A. For i = 1,2, let X; be the n x n matrix with a 1 in the i*® 
row and n‘* column. We consider four cases. 

Case 1: a(V) = 0 and a(M) = 0. M inherits laws i and ii from U and V. If 
a(U) = 0 then M also abides to laws iii and iv, and therefore is a uniquely pressable. 
This is contrary to assumption. We proceed under the assumption that a(U) = 1 and 
M is in violation of either law ii or iv. Since a(U) = 1 we have that the penultimate 
column of V is all ones above the diagonal. By law 727 the last column of V has at 
least n — 2 ones as well. Thus MW has exactly n—2 ones. Observe that if n is odd then 
V would violate law iv, therefore n is even. Finally, observe that M+ X satisfies all 


the uniquely pressable laws. 
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Case 2: a(V) = 1 and a(M) = 1. By law ii, since a(V) = 1, V contains all ones 
in its final column. Observe that M then inherits every law (i, i2, ii2,iv) from U. M 
is not a MNUP(n). 

Case 3: a(V) = 1 and a(M) = 0. Observe that the final column of M must 
contain n — 1 ones (by law ii on V). Therefore M inherits laws i and i from U 
and V. If n is even then law iv is violated, but M+ X, does not violate any laws. 
Otherwise, if n is odd then the only way that M can break laws 777 and iv is if the 
ante-penultimate column of M contains n— 2 ones. Observe that in this case M+ X, 
does not violate any laws. 

Case 4: a(V) = 0 and a(M) = 1. M violates law ii. Assume, by way of 
contradiction, that both M+ X, and M + X2 violate some law. Since M + Xo 
violates then it must be the case that the third row, n** column entry of M+ X, is 
a zero (otherwise it is the same argument as in case 2). Then V contains at most 
n — 3 ones in its final column. This implies, by law 777, that the penultimate column 
of V contains at most n — 3 ones and the ante-penultimate column contains at most 
max(2,n—5) ones. By law ii, the final column of U contains at most n — 3 ones and 
the penultimate column of U contains at most max(2,n—5) ones. These measurement 
extend to M and therefore M satisfies all four laws, a contradiction. 


Thus, we can make the following observation: 


Observation 6.9. If MM is a minimally non-uniquely pressable matrix then M + X, 


or M + X2 is uniquely pressable. 


Corollary 6.10. There are at most 3—(—1)"-3"'/4-! minimally non-uniquely press- 


able matrices on n vertices. 


Conjecture 6.11. Jf G is a minimally non-uniquely pressable graph, then every 
successful pressing sequence can be reached from another successful pressing sequence 


by a sequence of edits that involve at most four deletions or insertions. 
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Proposition 6.12. The Weak Pressing Game Conjecture on MNUPs is equivalent 


to the Weak Pressing Game Conjecture. 


Proof. Observe that if the Weak Pressing Game Conjecture is true, then it must also 
be true on MNUPs. We proceed to show that it suffices to verify the conjecture on 
MNUPs. Assume, by way of contradiction, that the Weak Pressing Game on MNUPs 
does not imply the Weak Pressing Game Conjecture. Let U be the instructional 
Cholesky matrix of a non-uniquely pressable graph G that is identity pressable but 
no other pressing sequence within four edits of the identity is successful. Let U be 
the instructional Cholesky of G. Observing that any 1 x 1 principal submatrix if U is 
uniquely pressable and G, by assumption, is not uniquely pressable, we may conclude 
that U contains a MNUP as a principal submatrix, choose one and call it M. Let 


A =UT'VU and consider the block construction of A: 


[ 4,| Ay | As | 


AZ| MT™M | Ay 


AT| AT | As 


Since the WPGC holds on MNUPs, there is a permutation matrix P such that 
P? M? MP is in pressing order and P can be converted into the identity permutation 
by removing and subsequently inserting at most four rows and columns (simultaneous 


operations on the same indexed row and column). However, 


T 
I} 0 {0 U, | Up | Us U, | Ug | U3 I|0}0 


0| P? | 0 0 |M|U, 0 | M|U, Oy. reese 
0} 0 [IL O> | 0b | Us 0 | 0 | Us 0; 0/2 
A, A2P A3 


PTA, | P?M?MP | P? Ag 


A3 AgP As 
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which can be successfully pressed by block-pressing A,, followed by block-pressing 


M?M (in the order given by P) and then block-pressing A;. However we now have a 


new pressing order for G that is within four edits of the identity, a contradiction. 


Conjecture 6.13. The MNUPs satisfy the Weak Pressing Game Conjecture. 


In particular, we conjecture that if M is a MNUP(n) it can be pressed by (exactly) 


one of the following sequences: 


e (n—1,n,2,3,4,...,2 —2,1,2), if M[1,n] = M[2,n] =0; 


(€,€+1,n,2...,€—-1,2+2,...,n—1,1), if 24 n is the first column other than 


1 with an odd number of ones; 


e (n,2,...,k —1,k+2,...,n,1,k,k +1) if M satisfies that M[i,n] = 0 if and 


only if 2<i<k-+1 for some even integer k; 


e (n,2,...,2—1,1) otherwise. 
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CHAPTER 7 


FUTURE DIRECTIONS AND OPEN QUESTIONS 


In this thesis we have explored several problems that are related to pressing sequences. 
Some of them we have been able to solve, many of them we have explored, and a few 
we have set on the back-burner. What follows are some future directions for this line 
of research as well as some open questions. We begin by addressing the elephant in 


the room. 


Question 7.1. Is the pressing game conjecture true? (Conjecture: Every successful 
pressing sequence can be reached from every other one by a sequence of edits that 


involve at most four deletions or insertions. ) 


A very small first step in the direction towards proving the pressing game conjec- 
ture would be to show that the Weak Pressing Game Conjecture holds up. Should 
the pressing game conjecture be false, it may still be true that successful pressing 


sequences can be reached by a relatively small edits. 


Question 7.2. Is it true that for a (full-rank) simple pseudo-graph on n vertices 
every successful pressing sequence can be reached from every other one by a sequence 
of edits that involve at most p*(n) deletions or insertions, where p*(n) is polynomial 


or less. 
Another very open and quite interesting (at least to me) question is: 


Question 7.3. What is the time complexity of determining the number of successful 


pressing sequences of a simple pseudo-graph? 
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We have shown that the complexity is no more than exponential. While this is an 
improvement over the previous upper-bound of super-exponential (derived from linear 
extensions), it is far from being computationally practical. Thus, if the complexity 


cannot be lowered to polynomial, we propose the following question: 


Question 7.4. Can the number of successful pressing sequences of a simple pseudo- 


graph be approximated near-uniformly in polynomial time? 


An affirmative to the previous question would likely lead to a near-uniform, effi- 
cient sampling algorithm for successful pressing sequences. 

In Chapter 2 we discussed enumeration of Cholesky roots for the zero matrix. As 
a corollary we were able to discuss the number of distinct Cholesky factorizations of a 
matrix (over F2) that is in leading principal minor non-negative form. This count only 


offers a lower bound for the number of distinct Cholesky factorizations of a general 


symmetric matrix over Fo. 


Question 7.5. How many Cholesky factorizations are there for a symmetric matrix 


over F,” Can we generalize the results for F,? 


In Chapter 3 we demonstrate that the successful pressing sequences of an OSP- 
graph are the linear extensions of a set of posets that arise from the instructional 
Cholesky roots of the graph. An autonomous graph has the property that its suc- 
cessful pressing sequences are all linear extensions of a single poset. In particular, 
in the autonomous case, this poset can be viewed as the intersection of all of the 
successful pressing sequences of the graph (interpreted as linear extensions). In the 
case that the OSP-graph is not autonomous then the posets are the intersections 
of pairwise disjoint families of successful pressing sequences. Thus, we have that if 
G is an OSP-graph then the instructional posets of G partition “(G) into disjoint 
sets $1, S2,..., 5, satisfying that LE(S;) = S; for each i € [k]. Observe that this 
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partition is not sufficient to determine the instructional posets of a graph since, for 


example LE (Q{o}) = {co}. 


Question 7.6. In general, how many distinct partitions of “(G) into disjoint sets 


{.5;}; exist such that LE ((|S;) = S; for each 7? 
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APPENDIX A 


SOURCE CODE AND EXAMPLES 


A.1 SAGEMATH CODE 


I have appended this section to preserve some of the SageMath code that we use 
to explore pressing sequences. Much of this code can be directly attributed to Josh 
Cooper, and the less elegant code is a result of collaborative work between Peter 
Gartland, Erin Hanna and I. Blakeley Payne (nae Hoffman) has also created some 
code to that is similar to what has appeared here - while it does not appear here, it 
is likely that her work inspired some of our work (thanks Blakeley!) 

The symbol “#” is used to denote a comment, and “<<<” is used following a 
line-break and means that the text which it precedes should be appended to the 
previous line. Throughout the code a finite field GF'(q) is assumed, prime power “q” 


should be defined prior to executing such codes. In the sections that follows I have 


included some sample outputs. 


HH 
# Input0: a vector ‘‘V’’ of length n(n+1)/2 
6 6 pe) 


# Inputl: the corresponding integer n 


# Output: an n x n upper—triangular matrix. 


Lip df pp ty py fy ty py Of pO Lip) J} fy Lip tpt py ty fp ty fp tp Of fp Lip ty fp py fy fp fy fy fp fy fy fit} fp ff ff 


def triangularize(v,n): 


loc = 0 
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outlist = [[0 for i in range(n)] for j in range(n)| 
for i in range(n): 
for j in range(i,n): 
outlist |i || 7) -= +] Loe] 
loc + 1 


return Matrix(GF(q),outlist ) 


# Input: a positive integer n 


# Output: 


list of upper—triangular matrices whose square is 


# the zero matrix 


Lip tp py pp ty fp py py OS Lt} Lt} 


Lift) ff py} ff py Lt} Lip d/) ff py pf fy ty ff fy Lift) ff fy} ff} 


def squarerootsofzero(n): 
size = int ((n+l)*n/2) 
“= matrix(GEeq) .n, ns 0) 
outlist =. [| 
for i in range(q size): 
S = triangularize(Integer(i). digits (base=q, 
<<< padto=size ) ,n) 
ibis 2. == Ls 
outlist .append(S) 


return outlist 


# Input: a positive integer n 
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# Output: list of upper—triangular matrices that when 


# multiplied by its transpose yields the zero matrix 


def Choleskyrootsofzero(n): 


size = int ((n+l)*n/2) 
A= matrix(GF(q)-;n3m.0) 
outlist =. |] 
for i in range(q size): 
S = triangularize(Integer(i). digits (base=q, 
<<< padto=size ) ,n) 
if S.transpose()*S = Z: 
outlist .append(S) 


return outlist 


# Input: a square matrix A 


# Output: True if the matrix is presented in leading 


# principal non—negative minors form, False otherwise. 


def is_LPN_k(A): 


n = len( list (A)) 
k = rank (A) 


for j in range(1,k+1): 


M = A. matrix_from_rows_and_columns(range(j),range(j )) 
if not M.is_invertible (): 


return False 
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for j in range(k+1,n+1): 


M = A. matrix_from_rows_and_columns(range(j),range(j)) 


if M.is_invertible (): 


return False 


return True 


# Input: an invertible square matrix A 


# Output: True if the matrix is presented in leading 


# principal non—negative minors form, False otherwise. 


def is LPN_full(A): 


n = len( list (A)) 
for j in range(1,n+1): 


M = A. matrix_from_rows_and_columns(range(j),range(j)) 


if not M.is_invertible (): 
return False 


return True 


Lip tp py ty fp ty py fp Sf pf Lf} Lip tf py py tf fp fy ff CY Lift} ff 


fit} ff fy ff {iL} Lip ft ff ff ff Lip ty) ff py} ff 


# Input: a matrix A satisfying is LPN_k(A)=True 


# Output: the ‘‘instructional’’ Cholesky root of A 


def iChol(A): 


C=copy (A) 
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n=len (list (C)) 
Bena trix (GF(2)4: my mg 0) ) 


for i in range(n): 


for k in range(i+1,n): 


C[k]= Cli, i]*C[i]*C[k, i]+C[k] 


return B 


# Input: A pair of posets, P and Q, 


over the same element set 
# Output: minimum edit distance between linear extensions 


# of P and Q 


def poset_edit_distance(P,Q): 


if set(P.list()) != set(Q. list ()): 
print ’The labels are not the same.’ 
return Infinity 
n = len(P. list ()) 
LP = map(lambda L: [P. list ().index(x)+1 for x in L], 
<<< P. linear_extensions()) 
LQ = map(lambda L: [P. list ().index(x)+1 for x in L], 
<<< Q. linear_ extensions ()) 
smallest_distance = Infinity 
for extP in LP: 
for extQ in LQ: 


perm = Permutation(extP ).inverse()* 
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<<<Permutation (extQ) 

d = n — perm. longest_increasing _ 
<<<subsequence_length () 

if d < smallest distance: 


smallest distance =d 


return smallest distance 


# Input: A symmetric n x n matrix A 


# Output: 


True if A has exactly one successful pressing 


# sequence, False otherwise 
Lp pp py pp ty py Ji dp ty ff fy 


{iL} Lf} {ft} ff Lt} Lt} Lift) ff py} ff py Lt} 


Lip d/) ff py pf fy ty ff fy Lip ty) ff fy} ff 


def is_unique_full(A): 
n = len(list (A)) 
presscount = 0 
for p in Permutations (n): 
PS pte imatrix) 
if is LPN_full(P. transpose ()*AxP): 
presscount += 1 
if presscount > 1: 


return False 


if presscount 


1: 
return True 


else: 


return False 


110 


Lip tp fy py py fy ty fy py fp 


tidy fp py tp fp fp Sf Lt} Lip ty py py pf fy fy ff fy 


Ld fp py fy py 
# Input:A symmetric n x n matrix A 


# Output: True if A has exactly two successful pressing 


# sequences, False otherwise 


def is_twice_pressable_ full (A): 


n = len( list (A)) 
presscount = 0 
for p in Permutations(n): 
P= pete mratrixt) 
if is LPN_full(P. transpose ()*AxP): 
presscount += 1 
if presscount > 2: 
return False 
if \ presscount-=—— 2: 
return True 
else: 


return False 


Lip ty py py ty fp ty fp ty Sf Sf) te dy dp py fp tf fp pf ff 


t/t} ff} LIP pe pe de ty ff fp fy LI fp fy fy py 
# Input: an invertible , symmetric n x n matrix A 


# Output: the number of successful pressing sequences of A 


def pressing number (A): 


n = len( list (A)) 
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presscount = 0 
for p in Permutations(n): 
PS p.to_ matrix () 
if is _LPN_full(P. transpose ()*AxP): 
presscount += 1 


return presscount 


# Input: a positive 


integer n and a positive integer k<n 
# Output: Permutations of {1,...,n} of length k extended to 


# length n by appending the missing elements in lex order 


Lip ty ff pf} pf fy Lift) ff fp ff tidy dy fy fp fy fy fy ff fp L}} 


Lip ty fp py fp fy ff fy Ly} Lip ty py pf fy fy ff fy Lift) ff py} fy 


def partialperm(n,k): 
permlist = Permutations (n,k) 
outlist = [] 
for p in permlist: 
plist = list (p) 
for j in range(1,n+1): 
if j not in plist: 
plist .append(j) 


outlist .append(Permutation( plist )) 


return outlist 


# Input: a symmetric n x n matrix A 
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# Output: the number of successful pressing sequences of A 


Lip ty fp py py ty ff fp fy 


Lett ft py tp fp fd pp tt) Sf Lt} Lip tpi py pf fy fy ff fy Lift) ff py} ff 


def pressingcount (A): 
n = len(list (A)) 
k = A.rank() 
presscount = 0 
for p in partialperm(n,k): 
Pop. to matric ()) 
if is_LPN_k(P.transpose()*AxP): 
presscount += 1 


return presscount 


# Input: A symmetric n x n matrix A 


# Output: the list of successful pressing sequences for A 


def pressinglistpartial (A): 


n = len( list (A)) 
k = A. rank () 
presslist = [] 
for p in partialperm(n,k): 
#print p 
P= pete, matrix) 
if is_LPN_k(P.transpose()*AxP): 
presslist .append(p) 


return presslist 
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# Input: An invertible symmetric, n x n matrix A 


# Output: the 


list of successful pressing sequences for A 


def pressinglist (A): 


len( list (A)) 

LIST =|] 

for p in Permutations(n): 

P= pete mratrixt) 

if is LPN_full(P. transpose ()*AxP): 


LIST. append(p) 
return LIST 


# Input: A graph G 


# Output: The adjacency matrix of G with entries in GF(2) 


def binaryadjmx(G): 


return Matrix (GF(2) ,G.adjacency_matrix ()) 


# Input: A matrix M and a non-negative 


integer m 
# Output: True 


if M has exactly m successful pressing 
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# sequences, False otherwise 


Lip tf fy py py fy ff fp py 


te dt ft py tp fp fp Sf Lt} 


eTNTETETETATSTATETSTRTATOTSTETST 
def has_number_of_PS(A,m): 
n = len(list (A)) 
k = A.rank() 
presscount = 0 
for p in partialperm(n,k): 
Pop. to matric ()) 
if is LPN_k(P.transpose()*AxP): 
presscount += 1 
if presscount > m: 
return False 
if presscount == m: 
return True 
else: 


return False 


# Input: a prime power q 


# Output: a list indexing the elements of GF(q) 


Lip dpe py py fp ty ty tt Sf 


Lip ty fp ty fp fp fy ft fp ff Lift) ff py} fy 


def FIELD(q): 


return GF(q). list () 
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# Input: positive integers n and r 
# Output: The number of symmetric matrices of rank r and 


# size n over GF(q) 


def NORMS(n,r): 


if min(n—1, r,n—r)<0: 
return O 


else: 


<<<for i in (0..r—1)) 


# Input: positive integers n and r and integer i in {0,1,2} 


# Output: the number of symmetric matrices of size n+l which 


# can be generated by increasing rank by i=0,1,2 


def NOM_gen(n,r,i): 


ib 0) 

return NORMS(n,r)*q r 
elif is=l: 

return NORMS(n,r)*(q—l)*q-r 
elif i==2: 


return NORMS(n,r)*(q~ (n+1)—q~ (r+1)) 
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else: 


return 0 


# Input: positive 


integers n and r 


# Output: DENOM(n,r) is used in the subroutines below 


def DENOM(n,r ): 


return NOM_gen(n—1,r,0)+NOM_gen(n—1,r—1,1) 
<<<+NOM_gen(n—1,r —2,2) 


# Input: positive 
# Output: 


integers n and r and integer 


probabilities used in subroutines below 


def prob_from(n,r,i): 


if i in range(0,3): 


return max(NOM_gen(n—1,r—i , i )/DENOM(n,r) ,0) 


else: 


return 0 


# Input: a symmetric matrix M 
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# Output: information about pivot columns of M used in 


# subroutines below 


def BASES(M): 


COMP=] 


n=len (list (M)) 


PIV=M. pivots () 


RPIV=M. pivot_rows () 
for i in range(n): 
F=matrix (GF(q) ,n,1,0) 
if i not in M. pivot_rows(): 
e=copy (E) 
e[i,0|=1 
COMP. append(e) 
cB=| 
for pan -PlY: 
CB.append( M.matrix_from_columns(|p]) _ ) 
P={ 
for p: in. PIV: 
P.append (p) 
return CB,COMP, P 


# Input: a symmetric matrix M 


# Output: a vector in the span of M 
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def inspanvector (M): 


L = BASES(M) [0] 


V = matrix(GF(q), len(list (M)),1,0) 
for i in range(len(list (L))): 
r = FIELD(q)[randint (0,q—1)] 
V=V+ raL[i] 


return V 


# Input: a symmetric matrix M 


# Output: a vector not in the span of M 


Lift) ff ty} ff fy Lip t/ ff ff ff Ld ty) fp ty ty fp ty fp ty Of fp Of Lt} Lip d/) ff py pf fy ty ff fy Lip f/f} py} fy} 


def notinspan(M): 
l 


len (BASES(M) [1] ) 
Y = matrix(GF(q), len(list (M)),1,0) 


s=ql-1 
it. «=i 

w= 1 
else: 


w = random_between(1,q 1-1) 
W= w.digits(q, padto=1) 
for i in range(1l): 
Y = Y + FIELD(q) [W[i]]*BASES(M) [1] [i] 


return Y 
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# Input: a symmetric matrix M 


# Output: a vector that is the sum of a randomly chosen 
# vector in the span of M and a randomly chosen vector not 


# in the span of M 


Lip d/ py py tf fp py fp ty Of fp Lip ty ff ty} pf py Led fp dy ty fp ty St OO 


def coolvector (M): 


return notinspan (M)+inspanvector (M) 


Lip ty ff py pf fy ff fp fy fe dy dt py fy fp yp Of Ly} Lip ty pe py pf fy py tO Sf ff 


# Input: a symmetric matrix M and an integer k in {0,1,2} 
# Output: a symmetric matrix of size one larger and rank 
# increased by k 


# Witty exceptions courtesy of Erin Hanna 


def rankinc(M,k): 


if. k==2) 
r = M.rank() 
if r = M.nrows(): 
print "I can’t increase the rank by 2!" 
else: 


Vv 


coolvector (M) 

p = M.rows() 

p[0:0] = (v. transpose ()) 
P= matrix(p) 


b = v.rows() 


F220 


b[0:0] = [(vector(GF(q), [9]))] 
C= matrix ('b) 

G = P. transpose () 

g = G. rows () 

g [0:0] = (c. transpose ()) 

OMG = matrix(g). transpose () 
BLAH = GF(q).random_element () 
OMG[0 ,0] = BLAH 

return OMG 


n = M.nrows() 
eta = [| 
for i in range(n): 
r = GF(q).random_element () 
eta.append(r) 
v = matrix(eta) * M 
p = M. rows () 
p [0:0] = (w) 
P.= Imatsix (p’) 
b = v. transpose ().rows() 
b[0:0] = [(vector(GF(q), [9]))] 
C= metre (b,) 
G = P. transpose () 
g = G. rows () 
g[0:0] = (c. transpose ()) 
OMG = matrix(g) 


NY = matrix(eta)*v. transpose () 


Lt 


iE: Kea 05 
OMG[0 ,0] = NY[0][0] 
elif k==1: 
x = NY[0] [0] 
while x — NY[0][0]: 
x = FIELD(q)[randint (0,q—1)] 
OMG[0 ,0] = x 
else: 
"you no do this" 


return OMG 


Lift} ff} Lip dy fp py fp fp ty ty ft tf Lt} Lt} Lift) ff py} pf py Lt} Lip dy ft py pf fy fy ff py {ft} pf L 


# Input: integers n,r,q 
# Output: a uniformly randomly selected symmetric matrix of 
a size n, rank r with entries in GF(q), or a witty 


+f exception message 


HH 
def nrmx(n,r,q): 
if min(n—1,r,n—r) <0: 
return "dummie" 
elit a4 
tt re U: 
Mematrix (GF(q) .{[0]]) 
return M 
elif r==1: 


Mematrix(GF(q) ,{ FIELD(q)[randint (1,q—1)]]) 
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return M 
else: 
return "exception in n==1 loop' 
else: 
P = |[prob- from(n 1,0). prob. from(n;r 1) 
<<<prob _from(n 92) | 


X = GeneralDiscreteDistribution (P) 


j = X.get_random_element () 
M=nrmx(n-—1,r—j ,q) 
M = rankinc(M,j) 


return M 
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A.2 CHOLESKY ROOTS OF ZERO 


2 x 2 square roots of zero: 


0 
0 
2 x 2 Cholesky roots of zero: 
0 
0 
3 X 3 square roots of zero: 
0 0 0 (ay 00 0 0 
00 0 4];,} 0 0 0],] 0 0 
0 0 0 0 0 0 0 0 
3 x 3 Cholesky roots of zero: 
0 0 0 0 1 0 0 0 
OO Oe is] Oe Oe ad Oe G 
0 0 0 0 0 0 0 0 
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4 x 4 square roots of zero: 


0 0 


0 


0 0 0 


0 


O° ol. 1.008 


0 
0 


0 
0 


0 
0 


0 
0 
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4 x 4 Cholesky roots of zero: 


0 0 


0 


0 0 0 


0 


Le 


0 
0 


1 
0 


0 
0 


0 
0 
0 
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A.3 LARGE UNIQUELY PRESSABLE MATRICES 


In Chapter 2 we discuss the set of uniquely pressable graphs by characterizing their 
Cholesky matrices. One can construct any uniquely pressable matrix by appending 
a column of full-weight to the end (and completing the matrix to be square), or by 
appending a new first row and writing a 1 in the entries of the first row exactly 
when the column previously had odd weight. We refer to to these processes as right- 
appending and left-appending (respectively). If n is even and U is a uniquely pressable 
nxn matrix, then the process of left-appending and then right-appending commutes. 
That is, all the (n+ 2) x (n+ 2) uniquely pressable matrices can be reached, without 
repetition, by performing the three following operations to the set of n x n uniquely 
pressable matrices: Left append twice, Right append twice, Left/Right append. Using 


this and the SageMath code below we generate random uniquely pressable matrices. 


def LeftAppend(1): 
w=[1] 
w. extend (1) 
for i in range(1,len(w)): 
if w[i]%2==1: 
w[ij=w[ij4l 


return w 


def RightAppend(1): 


x=copy (1) 
x.append(len(1)+1) 


return x 
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def RightRight (1): 
return RightAppend(RightAppend(1)) 


def LeftLeft (1): 


return LeftAppend(LeftAppend(1)) 


def RightLeft(1): 
return RightAppend(LeftAppend(1)) 


def LeftRight (1): 
return LeftAppend(RightAppend (1 )) 


def RandomUniqueEven(n): 
it a=): 
eutpub S| 1.2, 
else: 
i=randint (1,3) 
ae ee 
output=Left Left (RandomUniqueEven (n—2)) 
élit: T= 2: 
output=Right Left (RandomUniqueEven (n—2)) 
elif i==3: 
output=Right Right (RandomUniqueEven (n—2)) 


return output 


def SendUniqueToMatrix(L): 
n=len (L) 
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M=matrix (GF(2) ,n,n,0) 
M\Ecopy (M) 
for i in range(n): 
for j in range(L[i]): 
MM[i-j ,iJ=1 


return MM 


def RandomUniqueCholesky (n): 
if n%2==0: 
return SendUniqueToMatrix (RandomUniqueEven (n) ) 
elif n%2==1: 
1=RandomUniqueEven (n—1) 
i=randint (1,2) 
Li teat 
l=LeftAppend (1) 
elif i==2: 
1=Right Append (1) 


return SendUniqueToMatrix (1) 
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Figure A.l1 Random Unique Cholesky and Adjacency Matrix. N = 20 


Figure A.2. Random Unique Cholesky and Adjacency Matrix. N = 100 
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Figure A.3. Random Unique Cholesky and Adjacency Matrix. N = 500 


0 200 400 600 800 1000 1200 1400 0 200 400 600 800 1000 1200 1400 
of: * : : : : : * of ; ; ; ; : : : 
200 200 + 
400 400 5 
600 7 600 4 
800 800 4 
1000 + 1000 + 
1200 4 1200 + : 
1400 + 1400 + 
———— a 


Figure A.4 Random Unique Cholesky and Adjacency Matrix. N = 1500 
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